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FY 1977 FOUNDATION EXPLORATORY DEVELOPMENT PROGRAM 



I. Introduction 



The principal thrust of the exploratory development program 
at the Naval Postgraduate School (NPS) stems from the School's 
mission : 



To conduct and direct advanced education of 
commissioned officers, and to provide such 
other technical and professional instruction 
as may be prescribed to meet the needs of the 
Naval Service; and in support of the foregoing 
to foster and encourage a program of research 
in order to sustain academic excellence. 

In fulfillment of the exploratory development program 
objectives and within the above constraints, the Naval Post- 
graduate School 

Conducts exploratory development (6.2) deriving 
from scientific program areas or in other areas 
specifically requested by the Navy. 

In addition, NPS performs exploratory development, where 
uniquely qualified, for other agencies of the Department of 
Defense and, in defense related efforts, for other Federal 
Government agencies. NPS also furnishes consulting services 
for the Navy and, where specifically qualified, for other 
agencies of the Department of Defense and in defense related 
efforts for other Government agencies. 

II . Background 



The Navy has developed the Naval Postgraduate School as an 
academic institution which uses university educational methods 
to address the special graduate education requirements of the 
Navy. The Superintendent of the Naval Postgraduate School is 
a Rear Admiral of the line of the Navy. He is supported by a 
dual management structure, part military and part civilian. 

The faculty, mostly civilians, is responsible for the academic 
programs and, in support of these, conducts an active research 
effort. The military staff specifies the educational needs of 
the Navy, in terms of curricula, and provides administrative 
logistic support. 



At the Naval Postgraduate School, as in other academic 
institutions, the faculty is organized into departments. The 
department represents a resource center of faculty members 
with allied disciplinary specialization. Currently, the 
departments at the Naval Postgraduate School include: Computer 
Science, Mathematics, Administrative Science, Operations 
Research, Physics and Chemistry, Electrical Engineering, 
Meteorology, Aeronautics, Oceanography, Mechanical Engineering, 
and National Security Affairs. Inter-disciplinary groups 
which have effective departmental status include Electronic 
Warfare, Anti-submarine Warfare, and Command, Control and 
Communications (C3) . 

Development projects are largely conducted by the individual 
faculty members on a project basis. Projects typically originate 
from proposals prepared either by individual faculty members or 
by groups of faculty members from the same or different depart- 
ments . 

The development program is divided administratively into 
two parts. First is the Sponsored Development Program . This 
program includes projects awarded by sponsoring agencies to an 
individual principal investigator. The principal source of 
funds is the various commands and laboratories of the Naval 
Material Command. The sponsored program constitutes about 80% 
of the total NPS development program. The Foundation Develop- 
ment Program is based on a grant from the Chief of Naval 
Development. This program is administered internally by a Researcl 
Council, made up of faculty members appointed by the Provost. 

The Dean of Research is the Chairman of the Research Council. 

The Council meets periodically to review faculty proposals, 
allocate funding, and review results of completed projects. 

Appendix I identifies each project by title and category 
or type funding (Exploratory Development is 6.2 research). 

The 6.2 category was funded by the Chief of Naval Development, 
Washington, DC, 20360. 

This report summarizes the FY 1977 and FY 197T Foundation 
Exploratory Development Program. 
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Title : 


Investigation of a Sequential Test Median 


Investigator : 


D. R. Barr, Professor of Operations Research 


Sponsor : 


Foundation Research Program (6.2) 


Objective : 


To assess properties of a sequential test based 
on sample medians, including expected sample 
size and operating characteristics. 


Summary : 


Empirical evidence was obtained through simu- 
lation which demonstrates potential usefulness 
of the proposed test. This evidence suggested 
the sample number has asymptotic geometric 
distribution; efforts were undertaken to 
demonstrate this analytically, but to date have 
not been successful. 


Publications 
Forthcoming : 


It is anticipated that a journal article will 
be published on this material as soon as one or 
two open questions about the procedure are 
settled . 


Theses Directed: 


C. Y. Ching, "An Investigation of the Sequential 
Median Test," M.S. Thesis, NPS (1978). 

S. H. Wie , "A Comparison of Two Sequential CEP 
Tests," M.S. Thesis, NPS (1978). 



FURTHER REMARKS ON THE SEQUENTIAL MEDIAN TEST 



1. INTRODUCTION 

Testing hardware systems to determine whether they satisfy 
operational specifications is frequently an expensive and 
time-consuming task. In some circumstances, it is possible 
to decrease the consumption of resources in testing by using 
sequential test procedures. A sequential test requires that 
limits be placed upon the desired precision of the test 
(probabilities of errors) , and that the alternatives under 
test (hypotheses) be clearly specified. Such a test with 
two hypotheses is conducted by taking one observation (test 
value) at a time and, after each observation, making one of 
the decisions 

a. Terminate and reject the null hypothesis; 

b. Terminate and reject the alternate hypothesis; 

c. Take another observation. 

Since a sequential procedure requires that a decision 
(a - c above) be made after each observation, it is necessary 
to be able to obtain the observation value before another 
test is run (assuming it does not terminate at that point) . 

In the case of testing a system where several test runs may 
be scheduled in a relatively short span of time, this means 
that the measure being used in the test must essentially be 
obtainable in real time, without appreciable processing 
delays. Similarly, some "field computation," again without 
long delays, may be required in order to reach one of the 
decisions a - c. 

2 . Description of the Sequential Median Test 

Let the components of an IID sequence {X. , X_ , . . . } be 
observed sequentially, and let M. denote the sample median 
after 2k - 1 observations have been made. We consider a 
sequential test of H^: X^Fq vs. H^: X ^ F^ £ Fq based on 
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the statistics M k ; k = 1, 2, 3,... and sequences { (a k , r k ) } 
of continuation intervals whose end points are "acceptance 
numbers" a k and "rejection numbers" r k . We assume Fq 
and F. are absolutely continuous. The basic idea is to deter- 
mine r k ) } such that, if the test terminates when M k 

is not between a k and r k (with appropriate terminal action) , 
the procedure will have prespecified size a and power 1 - 6. 

The test we consider is based on medians of odd sized 
samples because the distributions involved are simpler, from 
a computational standpoint, than those for arbitrary sized 
samples. Even so, the situation is challenging, primarily 
because the sequentially observed sample medians are not 
independent. In previous work, we have developed a numerical 
method of determining the acceptance and rejection numbers, 
and have applied the method to a CEP problem like that intro- 
duced above . 

3 . An Application 

An example where such a procedure is useful is in testing 
air-to-ground weapon systems. A standard measure of perfor- 
mance for such systems is the "circular error probable" (CEP) , 
which is defined to be the median radial miss distance. 
Frequently, it has been assumed that the weapon (say, bomb) 
impacts under fixed drop conditions follow a bivariate normal 
distribution, and that by rotating from the "target plane" 
to a "normal plane" (perpendicular to the bomb trajectory at 
impact) , the distribution is approximately circular normal 
(bivariate normal with mean 0_ and covariance matrix a 2 1) . 

Under these conditions, the radial miss distances are Rayleigh 
distributed, and the familiar formula CEP =/TZn2 o 1.1774 a 
holds. There has been increasing citicism in recent years 
of using radial miss distributions based on bivariate normal 
impact distributions in the underlying plane. Common criticisms 
cite poor fit of the resulting radial model in the tails, 
especially in the upper tail due to a relatively high frequency 
of large miss distances in many data sets . Since there is 
apparently a serious question about the validity of such 
models in a substantial number of cases, tests of CEP by 
corresponding traditional tests of a 2 from an underlying 
normal model may be of doubtful validity. On the other hand, 
wide use continues to be made of models based on bivariate 
normal impact distributions, with corresponding faith in the 



implications of these models. The sample median, therefore, 
seems to be an attractive alternative in such cases, because 
it is generally more robust with respect to poor fit in the 
tails than are statistics sufficient for a 2 . In addition 
to its robustness, the proposed procedure has the feature 
that it can be used with truncated samples with little or no 
modification, unlike some of its competitors. The procedure 
is flexible, and can be adapted to meet various constraints. 
For example, one can accommodate taking observations in 
batches, also the procedure can be modified so as to preclude 
termination before some specified number of observations 
have been taken. 

The procedure is potentially of great usefulness for 
application to problems within DOD . For example, it should 
prove useful in connection with test and evaluation of a 
variety of weapons systems, as well as equipment components 
testing and quality control applications . 

4. Results 



Progress was made on the problems of determining asymptot 
properties of this procedure; even so, many challenging prob- 
lems remain. - The author continues to devote effort to this 
problem area as time to do so becomes available. 



Title: 



Investigator: 

Sponsor: 

Ob j ective : 



Summary: 



Conference 
Presentations : 

Publications : 



Aerodynamic Stabilization of Gaseous Dis- 
charges 

Oscar Biblarz, Associate Professor of Aero- 
nautics 

Foundation Research Program (6.2) 

The main objective is to define practical 
aerodynamic means for stabilizing discharges 
for electrical laser. Other applications in 
cold plasma chemistry are also envisaged. 

A particular objective was to become ac- 
quainted with supersonic flow and vortex 
flows as stabilizing agents and, furthermore, 
to examine the theoretical aspects of aero- 
dynamic stabilization. This is one part of 
a continuing project. 

The work undertaken complements the activi- 
ties in turbulence stabilization which have 
been underway at the Naval Postgraduate 
School. First-hand experience with super- 
sonic flow and vortex-flow equipment was ob- 
tained. The problem of analyzing the effects 
of turbulence in a discharge was given some 
perspective, and the classical or Schottky 
solution for the positive column in a flow 
discharge was re-examined in the light of the 
stability problem. The program will continue 
with more studies with turbulence. 



0. Biblarz, E. Wassertrom and Y. Crispin, 
"Stabilization of Flow Discharges by Super- 
sonic Flows", presentation at the 30th Annual 
Gaseous Electronics Conference, 18-21 Oct. 
1977, Palo Alto, CA, (abstract in conference 
record and in bulletin of the American Phys- 
ical Society. 

O. Biblarz, J. L. Barto, and H. A. Post, "Gas 
Dynamic Effects on Diffuse Electrical Dis- 
charges in Air", Israel Journal of Technology, 
Vol. 15, No. 1-2. 

O. Biblarz, J. L. Barto, and H. A. Post, 

"Gas Dynamic Effects of Diffuse Electrical 
Discharges in Air", 19th Israel Annual Con- 
ference on Aviation and Astronautics, Tel 
Aviv-Haifa, March 2-3, 1977, (full paper on 



conference record) . 



Y. Khait and 0. Biblarz, "Influence of Turbu- 
lence on diffuse Electrical Gas Discharges 
Under Moderate Pressures". Paper submitted 
to the Journal of Applied Physics . 

0. Biblarz, "Aerodynamic Stabilization of 
Gaseous Discharges", Technical Report, NPS- 
67Bi77111, 1 November 1977. 



Conference presentation at 30th Annual Gaseous Electronic 
Conference, 18-21 October 1977, Palo Alto, CA 
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Stabilization of Glow Discharges by Supersonic 
Flows* . E. WASSERSTROM, Y. CRISPIN and 0. BIBLARZ, 
Technion, Haifa, Israel. — It is known that there is a 
strong coupling between electrical discharges and gas 
flows. In our work, this interaction is studied for 
supersonic flow 1 , where the gas temperature is strongly 
coupled to the kinetic energy of the flow. Experimental 
results of supersonic flow discharges in cylindrical and 
conical tubes are described. It is found that the flow 
stabilizes the discharge and reduces the gas temperature. 
The results of some numerical calculations of a simpli- 
fied model indicate that the computed exit temerature 
agrees rather well with our experiments. Our work is 
relevant to the stabilization of electric discharges in 
lasers . 



1. E. Wasserstrom, Y. Crispin, J. Rom, and J. Shwartz: 
"The Interaction between Electrical Discharges and 
Gas Flow," to be published in the Journal of 
Applied Physics. 



*Submitted by Oscar Biblarz 













INFLUENCE OF TURBULENCE ON DIFFUSE ELECTRICAL 



GAS DISCHARGES UNDER MODERATE PRESSURES 

Y. KHAIT and 0. BIBLARZ* 

Technion - Israel Institute of Technology 
Haifa, Israel 

ABSTRACT 

The problem of a corona discharge under the influence of a 
turbulent flow is considered. This problem represents an extension of 
the already complicated situation of discharges in very non-uniform fields 
in non-moving gases, and our approach is restricted to providing a quali- 
tative description of important phenomena. The geometry studied is that 
of pin anodes and a wire mesh (plate) cathode immersed in a flow where 
controlled turbulence is generated by the presence of grids and screens. 

The action of these is discussed as to their desirable effects on the 
discharge. The applicability of the theory of locally developed turbulence 
is mentioned alone with semi-quantitative effects of turbulence on the 
transport coefficients. The fact that stationary distributions cannot be 
expected to occur under our conditions points out the usefulness of a 
pulsed-discharge analogy. Estimates of the relative power consumption, 
turbulent-to-laminar , can be made with the information given herein. 



Permanent Address : 



Naval Postgraduate School, Monterey, California. 
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PROOF OF ARTICLE TO BE PUBLISHED IN THE 
ISRAEL JOURNAL OF TECH., VOL 15, No. 1-2. 



Gas Dynamic Effects on Diffuse Electrical Discharges in Air 



O. Biblarz. J. L. Barto and H. A. Post 
Naval Postgraduate School, 
Monterey, CA 93940 USA 



ABSTRACT 

Effects of velocity from flows with and without grid-generated turbulence on the power input to a diffuse discharge 
are presented. Density dependences and some electrodes are also reported. Air has been used because, 
convenience aside, it presents some of the chemical complexities of molecular laser mixtures. Improvement of the 
diffuse discharge power is noted as a result of turbulent flows at subatmospheric density conditions — as much as 
500 times the maximum no-flow power. Density changes were obtained with a diffuser; the discharge geometry 
ordinary mode consists of an upstream pin-rack electrode and a downstream wire-mesh electrode, typically 2 V cm 
apart. A model lor the interaction of the ionization/recombination region at the anode pin tips is presented which 
desmbes some of The observations of the behavior of breakdown voltage with the flow. The voltage just prior to 
the glow collapse decreases somewhat in the ordinary (or anode-to-cathode) flow mode and increases in the res erse 
and cross-flow modes. The current just prior to lhe glow collapse increases marked!) under the action of a 



turbulent flow. 

NOMENCLATURE 

,. e — diffusion coefficient for ions, electrons 

— ambipolar diffusion coefficient (lami- 

nar or turbulent) 

— physical cap 

— electrical held 

s — electric field ai breakdown (just prior 

to the glow collapse) 

— electric charge 

— current at breakdown (just prior to 

the glow collapse) 

c — current density for ions, electrons 

— characteristic length for axial changes 

— characteristic length for radial 

changes 

— number density for ions, electrons 

n * — charge density 

INTRODUCTION 

as dynamic stabilization has been demonstrated to 
: an important means of increasing the power into 
fFuse discharges in molecular gases (Shwartz and 
asserstrom, 1975; Brown and Davies, 1972; Bib- 
rz and Nelson. ’97d; Wasserstrom et a!., 19 7 6). 



n 


— net rate of charge production 


R 


— “hyperboloidal” pin radius 


T 


— radial coordinate 


u 


— convective field velocity • 


V 


— discharge voltage 


v b , V*. 


— breakdown voltage (just prior to the 




glow collapse), reference value 


v t 


— inception voltage 


X 


— axial coordinate 


z 


— ionization coefficient 


a 


— recombination coefficient (tw-o or 




three body) 


p, Pu 


— density, atmospheric value 


& 


— anode extension for no flow 




— mobility for ions, electrons 


T 


— characteristic time 



These gases are found in molecular lasers. MHD 
devices, atmospheric discharges, and others. Turbul- 
ence, vorticity, and supersonic flow with shocks ha\e 
been demonstrated to accomplish more than just 
bulk convective cooling of the discharge; an interac- 
tion between 'he flow held and the ions takes place 



u 



which somewhat affects the fundamental nature of 
the glow discharge. For example, the so-called 
corona inception voltage and the breakdown voltage 
of the glow are modified by the flow, while at the 
same time the current becomes considerably more 
homogeneous and greatly enhanced. 

Our work has focused on the effects of turbulence, 
an effort which involves several related questions. 
We have been searching for an eddy size distribution 
or spectrum of turbulence which most favorably acts 
on the discharge (Biblarz and Nelson 1974; - W^w - 

1974), i.e., for some 
optimum which represents the most effective coupl- 
ing of flow power to discharge power. We have 
found eddy frequencies between 10 : and 10' Hz very 
favorable. The convective velocity itself has a 
stabilizing effect on the discharge, the understanding 
of which turns out to be important in the sorting out 
of the action of turbulence. The observed increase in 
current and power with velocity is nearly propor- 
tional to the square of the convective velocity 
( Nygaard , 1965; Biblarz , 1974). Colateral to these 
studies has been our interest in the effects of density, 
interelectrode spacing, and electrode orientation. 
Since our experiments arc carried out at pressures 
higher than optimum for discharges desirable for 
laser applications, it is significant to know the effects 
of reducing the density together w'ith changes of 
intei electrode spacing. The orientation (and to some 
extent the polarity) of the pin-to-plate discharge 
geometry with respect to the flow field has a bearing 
on discharge activity, and such knowledge turns out 
to be important in the understanding of gas dynamic 
effects, as will be shown in this paper. Finally, our 
experiments have been conducted with atmospheri- 
cal air where chemical complexities in the discharge 
are inherent; since we are observing relative changes 
in discharge performance and since air itself is a 
difficult discharge medium, we anticipe our results to 
be relevant to molecular lasers. Oxygen, for exam- 
ple, is used to increase laser power in CO-N 2 laser 
mixtures (Keren el al . , 1975); also, w'e expect our 
work to show qualitative trends in CO-mixtures. 

The plane in the positive point-to-point geometry' 
is represented in our experiments by a w'ire mesh 
which minimizes blockage of the flow'. Even though 
an array of pins is used to fill the interaction volume. 



the anode is analyzed on the basis of the activity at a 
single pin with some statistical allowance for the 
behavior of the set ( Bario , 197b). This type of 
discharge is “anode controlled” since the cathode 
plays only a secondary role ( Nasser , 1971). When the 
electric field has reached a sufficiently high value, an 
ionization region develops in the neighborhood o; 
the pin tips and it is this region that has the most 
profound effect on the discharge characteristics. The 
flow interacts with this region and a geometrical 
model which describes the anode extension, and 
which is discussed in this paper, seems to reflect 
some of the observed features of the discharge. 

The flow can interact with discharge phenomena 
that occur in comparable times. Table I indicates a 
few relevant characteristic times. Some values are 
obtained using a characteristic length of one mm, 
which is the order of magnitude of the observed 
anode region. 



TABLE I 

CHARACTERISTIC TIMES 





Characteristic 


Time (s) Reference 


Dilfusion (laminar) 


\0 l 




Dtlfusion (turbulent) 


10“ 4 




Flow time (1 0-1 00 m/s) 


10- -10- 




Turbulent eddies 


10 1 -10 5 


Biblarz and Nelson, 197 


Current pulsing 


5x 10 4 - 2x to 


* Nasser , t971 


Enorg) transfer, 






electrons to atoms 


to- -itr* • 


Meek and Craggs, t953 


Recombination 


10 ' 5 -10 * 


Meek and Craggs, 1953 


Maximum positive ion 






drift (£\= Itr V/m) 


10 3 -10 * 




Negative ion 






formation (CL) 


10 * 


Brown, !96b 


Streamer propagation 


!0- 


Nasser , !97t 


Maxwellian distribution 


to- -to "* 


.Meek and Craggs , 1953 


Excuaiion 


io-- 


Meek and Craggs, t953 


tom/ation 


10 


Meek and Craggs , 1953 



As may be evident from Table f, gas dynamic 
phenomena have times comparable to the recombi- 
nation time, current pulsing, and ion drift among 
others. The effect of flow on recombination itself 
may be overshadowed by the effect of flow on the 
current, as discussed in a later section of this paper. 
It is interesting to note here that current pulsing 



(Nasser, 1971; Loch . 1965), a statistical fluctuation of 
the current which manifests itself only at the onset' of 
the conyia and just prior to breakdown, has charac- 
teristic times close to those of the high end of the 
turbulence spectrum. Since the maximum positive 
ion drift time is not much faster than flow times, it 
may be surmised that within the anode extension, 
and particularly in regions of low electric field 
charge, convection should be significant. Diffusion 
times arc substantially lower than flow times and it is 
this difference which signals the possibility of using 
different characteristic lengths for diffusion and for 
convection in the equations (see Model and Addi- 
tional Results section). An important feature of 
discharges in air results from the dominant electron- 
attaching property of oxygen w'hich changes the 
electron energy distribution and electron density 
(Keren et a!. t 1975) and which cannot be affected by 
turbulence (see Table I). 

EXPERIMENTAL RESULTS 

Facility 

The experimental facility (Post, 1976) utilized con- 
sists of a plexiglass tunnel which is fed from an air 
compressor. Speeds in excess of 180 m/s can be 
obtained at relatively large flow rates with this 
facility. For subatmospheric operation, a diffuser is 
mounted after the test section. This diffuser was also 
made of plexiglass with an overall length of 1.04 m 
and wall divergence angle of 2.5 deg. No flow 
separation was observed in this diffuser. The test 
section houses the electrode system shown in Fig. 1, 
and the turbulence screen which is also evident in 
the photograph (Plate III attached to the pin rack). 
The intcrelectrodc spacing is adjustable but most of 
the result* reported here correspond to a spacing of 
2.9 cm. The anode consists of 39 pins forming the 
corners of rectangles 1.27 cm by 0.64 cm; the pin 
length is 1.57 cm. The test section cross-sectional 
area is 64 cm 2 . 

Most of the experiments were done with the 
electrode orientation indicated in Fig. 2a (ordinary 
mode). Some reverse flow and cross flow experi- 
ments were also carried out with the configuration 
shown in Figs. 2b and 2c. Although data^T" Loeb 
(1965) indicate that under no-flow conditions in air 
an optimum of the power into the discharge just 




prior to the collapse of the flow is obtained with a 
pressure of 500 torr (for a 2.9 cm gap), our systems 
w'as only capable of operating down to about 640 
torr (or p/p u or 0.85) because of the high blockage of 
the plate/electrode unit and because of chocking 
upstream of the electrodes. 

Electrical and fluid dynamic measurements were 
taken with conventional equipment (Biblarz and 
Nelson , 1974; Post , 1976). All measuring equipment 
was calibrated before use and corrections for humid- 
ity changes were made on the measured data (Post, 
1976), The screen used (Plate III) provided the best 
low frequency, high intensity eddies (Aunchrnan, 
1974; Post , 1976). Plate III consists of an aluminum 
plate laminated on phenolic with mixed arrays of 
holes of 1.0 cm and 0.3 cm diameter. The plate is 
mounted directly upstream of the pin rack and its 
area is smaller than the test section area to allow for 
a boundary layer “bleed”. In general we found that 
the pressure drop across ttie plate was a good 
indicator of its ability to stabilize the discharge 
( Aunchman , 1974). Turbulence intensities of the 
order of 10% are generated by Plate III. 






GRID 




CATHODE 



a) ORDINARY MODE 



U 



CATHODE — 



ANODE 



b) REVERSE FLOW 



GRID 



- > ' 



CATHODE 



ANODE 



c) CROSS FLOW 

Fig. 2 . Various orientations of electrodes and turbulence grids 
tested. 



Results 

Curves of power divided by the mass flow rate vs. 
the flow rate are indicative of flow effects (Shwartz 
and Lavte. 1975); any departure from the appear- 
ance of a rectangular hyperbola by the curves is to be 
attributed to flow effects. Figure 3 shows our experi- 
mental results for the ordinary mode in these coordi- 
nates; the power shown is the maximum diffuse or 
pre-arcing power. For the curves where the flow was 
turbulent the power was eventually limited by the 
power supply, whereas for the laminar flow curves 
the results were mass flow rate limited. The com- 
bined effect of reduced pressure and turbulent flow 
at maximum observed performance (Post, 1976) re- 
sulted in an increase of discharge powder of 75% over 
ambient laminar flow at the same flow rate. The 
initial downward trend of the curves in Fig. 3 indi- 
cates a constant power input; the slope then be- 




M ASS FLOW RATE THROUGH 
DISCHARGE (kg/5) 

Fig. 3. Power per mass flow rule vs. mass How rate through the 
discharge. 

comes positive as the flow interaction becomes sig- 
nificant (beyond 60 m/s). Pure convection seems to 
saturate in its ability to affect the discharge whereas 
turbulence does not. Error bars are shown in Fig. 3 
indicating the reliability of the measurements. 

Figure 4 displays breakdown voltage data as a 
function of velocity. Clearly, the major parameters 
affecting the breakdown voltage are the density and 
the velocity (fixed spacing); turbulence appears to 
have only a secondary effect. The decrease of V h 
with velocity is discussed in the next section in terms 
of the anode extension. Figure 5 is the corresponding 
plot for the breakdown or maximum diffuse current 
vs. averagCflow velocity; here, turbulence and veloc- 
ity appear to be the controlling variables with the 
density playing only a secondary role. The behavior 
of the maximum power follows closely that of the 
current, since the magnitude of the change in voltage 
is relatively very small. In addition to spreading the 
discharge at each pin (see Fig. 6), turbulence causes 
all of the pins to “light-up” more or less equally; 
under non-turbulent conditions the pins at the out- 
er periphery appear to be more active for the 2.9 cm 
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Fig. 6. Visual character of various glow modes without diffuser 



spacing; for spacings under 2 cm all the pins seem to 
be active under all conditions. Figure 7 shows photo- 
graphs of the various discharges. The no-flow photo- 
graph was a ten-second exposure. Here, the general 
appearance of the glow is that of a low- intensit) 
discharge extending between each pin and the 
cathode. The bright region at the bottom in the 
no-flow photograph was caused by a breakdown 
spark which took place during the time of exposure. 
In contrast, the laminar and turbulent flow photo- 
graphs were one-second exposures. In the laminar 
flow photograph, the individual pm discharges ap- 
pear to be more diffuse than in the no-flow case and 
an intense short narrow glow extends from the 
anode tips to the cathode. As evident from the 
turbulent flow photograph, the discharge here is 
even more diffuse with a wider bright region emerg- 
ing from the anode region; turbulent mixing seems 
to spread out the current and homogenize the cur- 
rent density throughout the interaction volume. The 
current-voltage characteristics of the diffuse dis- 
charge display the effects of convection and densitv 
reduction as typified in the information plotted in 
Fig. 8. 

Figures 9-12 will be discussed as part of the next 
section in which a model for the anode extension is 
presented. 




Fig. 7. Discharge appeal ance. 




Fig. 8. T>pical observed I-V shift due to diffuser induced den* 
sily reduction. 

MODEL AND ADDITIONAL RESULTS 

Model 

The observations depicted in Fig. 6 indicate de- 
finite affects of the flow on the radiative recombina- 
tion pattern at the anode tips. -We propose to discuss 
here an elementary model which is based on a 
geometrical interpretation of changes attributable to 
diffusion and convection within the anode region. In 
this region, which exists already without any flow, 



intense ionization and recombination take place. 
Even though the anode region includes the sheath, 
we shall assume the region to be quasineutral for the 
simple reason that the sheath is considerably smaller 
(Biblarz et al. y 1975) and least affected by convection 
due to its proximity to the wall (Dotson and Biblarz , 
1976). If there are changes in the anode region 
attributable to the flow, then we can think of them as 
changes in the effective physical gap which in turn 
will affect breakdown-gap related parameters such 
as V c and V b . In particular, it has been observed that 
the flow extends the anode region in the ordinary 
mode of operation and that this results in a decrease 
in the breakdown voltage. 

From Figs. 6 and 7 it is seen that the diffusion~field 
for the no-flow case is more or less spherical about 
the pin tip. Convection stretches out this diffusion 
region into an oblate spheroid and turbulence fat- 
tens the region considerably. 

Since we are disregarding effects of convection on 
the sheath, the anode extension can be described 
with the quasineutral equations ( Shwartz and Was - 
serstrom , 1975) where diffusion is ambipolar and 
where the effect of convection is included (Biblarz 
and Dolson , 1974; Sutton and Sherman , 1965). The 
species current equations are 



j, = en.ji.E + en U - eD.Vn,, 
j K = ~ en c fi c E > en e O - eD e Vn ei 



( 1 ) 

& 



here the symbols are defined in the Nomenclature, 
he total current is 7 = In the ambipolar 

igion, n, = n c - n and the basic equations that 
escribe the system become ^ 



V • E = 0, (3) 

V-h. c =eru (4) 



elutions to Eq. (3) would have to be ’‘patched” to 
listing sheaths and subsequently to the external 
rcuit. Equation (4) assumes steady flow; we shall 
eat turbulent flows as pseudosteady where the 
fleets of fluctuations appear in the form of greatly 
nhanced transport coefficients. The form of the net 
roduction term, ri, that is ususai in positive column 
nalyses ( Cobine , 1958; Jaeger el al ., 1976) includes 
mization by electron collision but two-body recom- 
ination (i.e., Zn - an 2 ); this is not appropriate for 
ear-atmospheric operation where three-body pro- 
ves dominate. We shall use, therefore, the three- 
ody form of Hinnov and Birschberg (1961) 



h - Zn - an 3 . 



( 5 ) 



he above set of equations is incomplete without a 
jitable energy equation ( Sutton and Sherman , 1965) 
hich yields the electron temperature and the gas 
imperature, if Joule heating is appreciable, and 
)me information regarding the flow field. If we 
liminate the electric field from Eqs. (1) and (2) in 
le usual way (Brown, 1966) we obtain 

n V • 0 + 0 • Vn - V • (DSn ) = Zn - an\ (6) 



'ere D a is the ambipolar diffusion coefficient. The 
hove equation can be specialized to a simplified 
ersion of the conditions for the ordinary mode (see 
ig. 2a). Assuming that U is only in the axial 
irection and constant and that diffusion is radial 
lainly and symmetric, we can write Eq. (6) as 




ld_ 
r dr 




Zn — an J . 



( 7 ) 



he above is a non-linear, second order paitial 
irterential equation which we shall nor attempt to 
alve here. It is useful, however, to suitably non- 
imensionalize (Kline, 1965) the equation so that 



characteristic parameters may be estimated. First of 
all, let us define n * such that 



Zn * — an * 3 = 0. (8) 

This corresponds to an equilibrium condition in 
which ionization and recombination locally balance 
each other. Note that n * is a function of the electron 
temperature through the Saha equation (a is also a 
function of electron temperature (see Hinnov and 
Hirschberg , 1961 or Mitchner and Kruger . 1973)). The 
parameter n * will be used for the characteristic 
density since n/n * ^ 1.0. As stated in the Introduc- 
tion, there are really two characteristic lengths in the 
problem, one for the x-variation which is primarily 
convective and another for the r-variation which is 
primarily diffusive. With this in mind, and introduc- 
ing D J as the ambipolar diffusion coefficient for the 
condition n* y we can rewrite Eq. (7) as follows: 



(JD 


i dn 


mu* 


\ZL.t 


' dx 


where 




n «-* n In*, 






x x / 






r «-» r/L n 



Dt*+DJD*. 

Since we have made the variables of order one, it is 
reasonable to assume that each group of coefficients 
that appears in Eq. (9) is also of order one. Hence, 
we can arrive at an estimate for the characteristic 
lengths: 



L =1^- 
* 2 ♦ 


(10) 


an 


L, = yj D If an * 2 . 


(in 



Within a cylinder characterized by a length L, and a 
radius L n we anticipate to find the anode region. The 
accuracy of the above estimates will depend on the 
extent to which the electron temperature is constant 
in the ambipolar region because we have inserted n * 
and D * through the differential. Note also that we 
have assumed a thermally stable system. 

Equations (10) and (11) are the desired interim 
i - ults. Equation (10) agrees with a solution of E'.q. 



(9) in which diffusion's completely neglected. Equa- 
tion (11) is consistent with work reported in Mitch- 
ner and Kruger (1973) where a solution of Eq. (9) 
without the convective terms is shown to be of the 
“boundary-layer type”. We can state now that the 
condition for which diffusion can be neglected in the 
■^-direction is one for which L, > Lr. Note that L , 
will be greater for turbulent flows than for laminar 
ones. 

It is clear that, other things being equal. L t will 
increase with increasing convection velocity. If 
this corresponds to a lengthening of the 
ionization/recombination region then our geometric 
model appears to be appropriate, fn particular, from 
Fig. 4, we have found that the breakdown voltage 
decreases with increasing velocity; Fig. 9 shows data 
additional to Fig. 3 for VJVw where is the 
no-flow value. In this figure, an eventual increase in 
\/ b is also apparent (a more appropriate variable for 
the ordinate is U / d because it collapses the data for 
various gaps into one curve). The eventual increase 
in breakdown voltage, if it really exists, is not easily 
explained. It is known that turbulence causes the 
electric field to jrise because of greater charge losses 
( Shwartz and Wasserstrom , 1975). This turbulence 
(either from the wake of the pins or grid generated) 
increases in intensity with increasing velocity, but it 
is not clear that results for a fixed column discharge 
apply to our case. Another possibility lies in t lie 
effects of convection on recombination; it is reasona- 
ble to expect that L x is shrinking because n* is 
increasing very rapidly; since the current in the 
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discharge is roughly proportional to U 1 and since 
appreciable electron heating may take place at the 
higher currents, we can perhaps attribute the rise in 
voltage to changes in recombination length arising 
from an inciease in the denominator in Eq. (10). 



Additional Results 

Figure 10 shows VJV data for the reverse and 
cross-flow modes. Here, according to our model, the 
anode extension which exists for no flow is eventu- 
ally, at high enough speeds, completely blown out of 
the interelectrode region leaving the physical gap as 
the smallest (or breakdown) gap. It is evident from 
the data in Fig. 10 that the voltage increase levels out 
at the high velocity end. These results are consistent 
with those of Biblarz and Nelson (1974) and Werner 
(1950). If d is the physical gap and 5 () the anode 
extension for no flow, then the maximum increase in 
the breakdown voltage is simply 
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Fig. Voltage ratio vs. U/U tn the ordinary flow mode. 



Fig. 10. Voltage ratio vs. velocity in reverse and cross flow 
modes. 



We are implying, of course, that the breakdown 
voltage is proportional to the distance between the 
anode or anode extension and cathode ( Barto , 1976; 
Nasser, 1971). The sheaths at these pressures are 
very thin. In order to arrive at an estimate of 8 0 we 
must note that it’should be a function of the electric 
field distribution outside of the anod^region. As- 
suming a hyperboloidal point-to-plane distribution 
(Barto, '1976) for the electric field in the neutral 
region we have arbitrarily chosen the anode region 
to lie where E ^ 20 kV/cm; this is a value of the 
electric field which is greater than the 15kV/cm 
needed for ionization by electron collision in air. 



Therefore 


. # 


e( ')-rL e - 


(13) 


where 




£ = 2V 
p R\og(4d/R) 


(14) 


and 5„ = x when E(x) = 20kV/cm. 





Using values calculated with Eq. (13) and taking 
the data for the leveling off value of VJV W from 
cross-flow and reverse flow measurements, we arrive 

■s 

at the information shown in Table II and in Fig. 11. 
It must be kept in mind that the data come from 
various velocities and levels of turbulence. In Eq. 
(14) we use R = 0.05 cm and V b = 7.8d + 4 (kV) as 
suggested by our experimental conditions (Barto, 
1976). The trend is well indicated by the calculations, 
and the agreement between theory and experiment 
appears to be reasonably good. This then is support- 
ing evidence for the anode extension model. 
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REVERSE ANDCROS* FLOW MAXIMUM VOLTAGE RATIOS 
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Fig. U. Summary' of maxium breakdown voltage ratio vs. gap 
length, reverse and cross flow only. 



Our observations during the cross-flow experi- 
ments revealed a host of interesting effects, and since 
these may have a bearing to £- beam lasers they will 
be briefly discussed here. Figure 12 is a summary 
schematic of our observations. The ions within the 
anode extension, 8, are affected by two orthogonal 
forces, namely, convective drag in the flow direction 
and the electric force in the cathode direction. The 
net effect is shown in Fig. 12 as a displacement to the 
right which acts to lengthen the interelectrode gap 
and to distort the local field at the anode tip region. 
The displacement of the anode region was such that 
apparently it could affect the downstream pins (1.27 
cm away!). Instabilities were observed in the second 
and third rows of pins (b and c in Fig. 12) in the form 
of minute but numerous transient arcs. These minute 
arcs were not strong enough to trip the power 
supply; they, however, spread into the first row of 
pins just prior to breakdown (when the power supply 
is actually tripped). Since breakdown streamers arc 
energetic enough to be uninfluenced by the flow, ihe 
displacement of the anode region shows up as an 
increase in the breakdown voltage, as discussed 
previously. Even while these numerous arcs were 
evident, a blue-violet glow region (region 1 in Fig. 
12) marked an upstream boundary for these arcs. 
This glow region was visible only from the first and 
second row of pins and it always “leaned into the 
wind”. Downstream from the last set of pins, a 
blue-green glow was evident the extent of which 
decreased with increasing velocity; *«■•*.•» ‘ n 

Fig. 12 was of the order of 5 cin. In B - row tr(19jtS) ii is 
shown to follow a recombination controlled pattern 
which is proportional to the inverse of the convec- 
tu* velocity. 



v V \ V •. 




Fig. 12. Schematic representation of cross-flow current distribu- 
tion observations. 



CONCLUSIONS 

The power into a diffuse discharge in air has been 
enhanced by 500 times by the combined action of 
flow, turbulence and reduced density. To some ex- 
tent, effects of velocity, turbulence, and density have 
been separated allowing individual scrutiny of each. 
The experiments have been undertaken in a rela- 
tively high flow facility with a “difficult” molecular 
gas mixture. 

By far the most pronounced effect is the enhance- 
ment in the magnitude of the current. Turbulence, 
for example, makes all the anode pins active, which 
is a difficult task to accomplish in practice with gaps 
over 2 cm (no individual ballasting of the pins was 
used); the increase of current, however, was more 
than the number of pins that were “activated”. 
Turbulence considerably enhances and homogenizes 
the anode region apparently allowing each pin to 
flow more substantial amounts of current. Of course, 
velocity also shows up as an important parameter, 
but there is no pure laminar flow at the pin tip 
because it is in a wake where vorticity is being shed; 
still, the current has been observed to possess a 
square or higher ( Barto , 1976) dependence on the 
convective velocity. It is interesting to note that 
density changes by themselves did not seem to alfect 
the results shown in Fig. 5; this is somewhat surpris- 
ing in the light of the dependence of the Townsend 
coefficient on the density ( Nasser , 1971; Cobine , 
1958). Figure o shows some density enhancement but 
not a substantial one. Turbulent flows have been 
seen to stabilize our discharges in which Joule heat- 
ing and thus thermal instabilities are not as substan- 
tial as in the laser discharges ( Drown and Davies, 
1972; Wdsscrstrom et al. y 1976); it may therefore he 
concluded that the elfect is a fundamental one. 



With the anode region model we have shown that 
it is possible to interpret the action of convection in 
terms of a physical extension or contraction of the 
anode region. This incidentally can also be argue to 
apply to results with parallel plate electrodes ( Batio , 
1976). In the ordinary flow mode, a minimum of 
VJVw is obtained of about 0.90-0.95 corresponding 
to a velocity of about 75 m/s; the subsequent in- 
crease in breakdown voltage is at the present not a 
conclusive result. The reverse and cross-flow modes 
saturate at a VJV ^ of about 1.25, which is post- 
ulated to correspond to the limit imposed by the 
physical gap. The results for the breakdown voltage 
indicate that a relation exists of the form 




where r is some relevant characteristic time. 
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I. INTRODUCTION 

In the high energy laser field, the electric discharge convection laser 
(EDCL) offers good efficiency as well as simplicity of handling relative to 
gasdynamic and chemical lasers. Progress towards higher laser outputs from 
EDCLs, however, has been impeded by the collapse of the glow discharge at the 
increased currents and pressures desired. This glow collapse (usually re- 
ferred to as arcing or constriction) represents the practical upper limit of 

the electrical power that may be pumped into the necessary non-equilibrium 

1-4 

processes. It has been shown that aerodynamic stabilization can be used 
very successfully in overcoming the glow-collapse limits. It should be 
pointed out here that the laser application is but one application of the 
relatively new interdisciplinary field of cold plasma chemistry. 

Aerodynamic stabilization is based on two separate effects usually re- 
ferred to as convection (laminar) and mixing. The purely convective effect 
is said to decrease the residence time of the particles through the discharge 
and thereby to counteract unstable mechanism that lead to the glow collapse^. 
The mixing effect, which may be due to vortex flows or to turbulence from 
boundary layers or grids, has been observed to have powerful stabilizing 
characteristics. Convection scales up with flow velocity whereas mixing 
depends on the intensity, spectral distribution, and location of the eddies 
or vortices . 
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II, SPECIFICS 



At the Technion, the author had the opportunity to participate in 

experimental work dealing with vortex-flow stabilization^ and with super- 

6 8 

sonic flow stabilization * . These techniques complement the work that has 
been underway at the Postgraduate School on turbulence stabilization of 
discharges. Some results of the Technion work have already appeared (see Refs. 
6-11) and more will be forthcoming. 

A significant portion of the author's time was dedicated to a study of 
theoretical aspects of discharge stabilization. The author was privileged to 
colaborate with Dr. Y. Khait in establishing a perspective regarding some of 
the author's previous results with turbulence stabilization. This effort 
resulted in the writing of a paper titled "Influence of Turbulence on Diffuse 
Electrical Gas Discharges Under Moderate Pressures," which is being submitted 
for publication. Other theoretical work involved a close look at the 
Schottky solution for the positive column in a glow discharge with the aim 
at establishing its use for a stability formulation. This effort was fruit- 
ful on two intermediate aspects of the problem and these are included here 
as Appendixes A and B. 

A proposal was written to the United States-Israel Binational Science 
Foundation^ which extensively outlines the progress of the entire effort. 

If awarded, the cooperation of the Technion and the Postgraduate School 
would be formalized (however, the budget would provide neither salaries for 
principal investigators nor equipment for the USA experiments) . 
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III. DISCUSSION 

The classical solution to the ambipolar diffusion of charges in the 
12 13 

positive column ’ of a glow discharge offers no hints as to how striations 
and constriction, among other important observations, manifest themselves in 
the discharge. As pressure and/or current are increased, instabilities grow 
in the discharge and it is important to have a sufficiently complete model 
of the physics which reflects these effects. In other words, before flow 
effects such as (laminar) convection and turbulence are introduced it is 
necessary to re-examine some assumptions inherent in the Schottky solution. 

A partial list of phenomena not reflected in the Schottky solution may 
be given as follows: 

1) effect of sheaths at the walls 

2) effect of electrodes and intermediate regions 

3) effect of recombination in the gas (pressure and composition) 

4) effects of the energy balance, including radiation and the elevation 
of the electron temperature 

The list can be rather long. We closely examined the energy balance because 
the instabilities are believed to be thermal disturbances (i.e., the tempera- 
ture becomes non-homo geneous) and recombination in the gas because it closely 
reflects the effect of pressure and composition. For simplicity, a two- 
dimensional (infinite parallel plate) geometry was considered, see Fig. 1. 

When E t is constant, V*J = 0 implies that J = J(y) or that 
n = n(y) because J^_ = 0 ; moreover, V x E = 0 implies that E = E(y) . 
The quasineutral equations for steady state may be written as^’^ 
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-Figure 1. Geometry of Positive Column 
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with boundary conditions 



at y = 0 



3n/3y = 3T/3y = 0 



at y = + a n = 0 and T = T 

' — r 



la the above formulation x is the axial coordinate and y the interplate 
coordinate, n the charge concentration in the ambipolar region, the 

electron temperature and T the gas/ ion temperature, N is the concentra- 
tion of the neutrals, the coefficient of ambipolar diffusion, K the 

thermal conductivity, the electron mobility, and a the ionization 

and recombination coefficients respectively, E the axial electric field, 
v ^ the collision frequency between electrons and neutrals, e the elec- 
tronic charge, p the pressure, and and m^ the mass of the electrons 

and neutral respectively. There are several assumptions inherent in this 
formulation, namely that N>>n , T^>T , and that p = constant . In addition, 

we are using a symmetry condition in the boundary conditions and we are 

12 13 

neglecting of wall space charges. Glow discharge data 9 , in addition, 
indicate that E^ as well as appear to be constant in the discharge at 

least before constricting. The electron energy equation, Eqn. 3, should have to 
reflect the energy dependance of the collision cross sections and, in addition. 
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non-elastic collisions — since we are interested in laser operation. 
However, we shall not deal with Eqns. (3) and (4). 

Using simple kinetic theory to describe the gas temperature dependance 
of the coefficients (D , p , K , etc.) and defining a stationary state 

cl 6 

where 

v^n* - an* 2 = 0 (5) 

or n* = v_^/a 

we can proceed to nondimens ionalize Eqns. (1) and (2). We shall, however, 
not use the plate spacing to render y dimensionless, but introduce two 
characteristic lengths, and . The gas temperature is made non- 

dimensional by the use of a wall temperature T q . The resulting equations 
are: 



, ,, /-,l/2 dn, 2 

d/dy (T ^) - n - n 


(6) 


IT D 

r - / e ao 

^ " V 2v io T 0 


(7) 


d/dy (T 3/2 - nT 


(8) 


/ T K 

T / 0 0 

T ~ Ven*u E t 2 

eo L 


(9) 



It seems reasonable to expect that in a conventional (unconstricted) 
glow, the two characteristic lengths be about equal to each other and that 

L n /a s 0(1) and L T /a « 0(1) (10) 
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This expectation is also required by the eigenvalue nature of Eqns. (6) and 
(8). It may be shown (see Appendix B) , that except for the equilibrium limit, 
Eqn. (6) is an eigenvalue equation, (the classical solution being a cosine 
when no recombination is present and when T « constant) . The nature of 
Eqn. (6) is discussed further at the end of this section. 

In the light of EqnS. (10) , it is interesting to look at some actual 
values of L q and . Particularly extensive information for helium is 
reported in Ref. 16, and this is used in estimating the values shown in 
Fig. 2, which are limited to the range of applicability of the equations 
in Ref. 16. An appropriate form for the coefficients is obtained by intro- 
ducing the gas density N q , 
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The above forms show "similarity" since LN is a function of E/N only. 

Examining Fig. 2, it is evident that L^N is not equal to L^N , certainly 

not for large E/N , wheareas L,^ follows closely data reported for helium 

12 13 

in the literature * (the fact that we are looking at parallel plates and 
the data is for tubes is immaterial) . Since in a discharge tube n and T 
share a common diameter, it is necessary to reconcile the discrepancy of the 
two curves. 

14 

It is possible, that the assumption which brings T g into Eqn. (1) 

and, therefore, into Eqn (11) is invalid. But, even if we were to eliminate 



7 



HELIUM 






ro 

i 




.. CD 





to 






ro 



CM 



i 



to 

I 



I 



I 



LOG (LN) 



39 



the factor T e / T 0 from the equation, the curves would diverge at the higher 
values of E/N . Constriction of the positive column is another possibility 
but in this case would be expected to be greater than for any given 

E/N and just the opposite seems to be true. Uncertainties in the estimate 
of a , the recombination coefficient, would affect the value of but 

this estimate is the one that follows nicely the experimental values. 

The method of reconciliation of the two lengths proposed here is based 
on the alternative that the equations previously written are incomplete and 
that, in particular, in Eqns. (1) and (7) should be replaced by a 

parameter (C^) which makes up for the difference in the two curves. This 
matter is discussed in Appendix A of this report. 

Consider now Eqn. (1) with constant coefficients. We proceed here to 
nondimens ionalize it in a conventional way 

- D n" ■ v. n - an^ (13) 

a i 

n’(o) = 0 and n(+ a) = 0 



Let n(o) = n and e = an /v. 

O O 1 

Hence we get 

- n" = n - en^ 

n* (o) = 0 and n(o) = 1 
n(+ a/ \j D /v.) = 0 

cl X 



(14) 



Now the Shottky solution is one for which e ■ 0 , namely. 



A »• ■** 



and n = n cos (•=■ y/a) 

O L 



(15) 

(16) 
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It can be seen therefore that the characteristic length L should be 



L = 




a 

tt/2 



(17) 



On the other hand, for the equilibrium limit e = 1 and the equation 

*2 



(18) 



is of the "Boundary Layer" type, i.e., it is not an eigenvalue problem, see 
Ref. 14 pp 149-51 for the equivalent problem with three-body recombination. 
The solution to Eqn. 18 in an implicit form is 



y = log 






| + (h + |) 1/2 



I - (h + \) in 



1.317 



(19) 



or explicitly 



3 

n = 2 n o 



1 _ e - (7 + 1.317) 



1 2 



1 + e 



- (y + 1.317) 



n 



( 20 ) 



which is indeed of the boundary-layer type. The intermediate cases, 

0 < e < 1 , are represented by elliptic functions, and this matter is dis- 
cussed in Appendix B of this report. 
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IV. CONCLUSIONS 

Some deficiencies of che Schottky solution have been re-examined with 
the intent of establishing a more useful solution for the description of 
the transition to the glow collapse. At least two major, recent attempts in 
the literature (Refs. 16 and 17) base the steady solution for the charge 
profile on the Schottky result and these attempts have not succeeded in pre- 
senting a complete formulation of the problem. The Technion approach^ is 
to base the perturbation on a physically complete solution of the ambipolar 
charge profile. This profile must include the ingredients of constriction 
and of the elevation of the electron temperature. 

As shown in Appendix A, the characteristic length scales for ambipolar 
diffusion and heat diffusion (conduction only) may be reconciled if in 

Eqn. (7) is replaced by where and C 2 is a function of 

E/N (i.e. , the dependance given by the difference of the two curves in Fig. 
2). What appeared to be a bonus, namely, an insight into the mathematical 
description of striations, needs further study in order for it to make more 
physical sense. The eigenvalue nature of the diffusion equation is not 
changed by the recombination term except at the limit when t - 1.0. 
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APPENDIX A - PROBLEM REFORMULATION 

Equation (1) in the text presupposes steady flow along with n = n(y) . 
Let us relax these suppositions and go back to the species flux equations 
which may be written as 

D 



I7 + V* [n(- u E - -r 2 - V n) ] = n 
<3 1 e n 



D. 

I 7 + V- [n(M. E - V n)] = h 
d t In 



(Al) 

(A2) 



Combining and eliminating the electric field we get the ambipolar diffusion 
equation in the usual manner which for constant coefficients is 



or 



3n/3t - D V n = h 
a 

3n/3t - D (£± + iV 
a W 8 y 2 / 



= n - an 



(A3) 

(A4) 



The above assumes two-body recombination as well. Now, using Eqn. (A3) to 
2 

eliminate V n from Eqn. (A2) yields 



an . "1 7 ' nE • 

3t 1 - D./D n 

1 a 



(A5) 



But since V*E = 0 in the ambipolar region 



3n i 

3t 1 - D./D 
x a 



E*V n = n 



(A 6 ) 



or 3n/3t + 



,1 - D i /D a ) f E L 3x + E y 3yy " V i n " an 



(A7) 
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Let us work with Eqns. (A4) and (A7) in linear form, i.e., neglecting 
recombination. If 

n(x,y,t) = f (y) g(x,t) (A8) 

Then, substituting into Eqn. (A4) 

f S t " D a ^ fg xx + gf,,) = V i fg ~ ••• (A9) 

2 2 

where f" = — | g^ - and g t = 3g/3t 



Dividing by fg both sides of (A9) 

1/g (g t - D a gxx ) - f f " = v. (A10) 

At a given E/N , is constant so we may presume that 

1/g (g_ “ D g w ) = constant (or - C ) (All) 

£ cL aJC X 

D 

and - ~ f" = constant (or C^) (A12) 

with ~ C 1 + C 2 = V i (All) 



The constants and C ^ are positive and defined so that 

as required by Eqn. (A12) , which now replaces Eqn. (6) since 

- D f" = f C. (A14) 

a l 



and L.. = \J D /C. 

f a l 
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Equation (A14) is the new Schottky equation. The constant appears in 

the companion partial differential equation 

*t ' D a 8** + S - 0 

Now look at Eq. (A7) . At the centerline 3n/3y = 0 and f(o) = f 

o 

so that 



at y = 0 



3n . __ i L 3n _ , 

3t + ll - D i /D a j' E L 3x i n ” "* 



Introducing n = fg we obtain 



g t f o + \1 - D./D / E L 8 x f o v i f o e 



i a/ 



(A16) 



(A17) 



We can cancel in the linear equation above and obtain an equation 

in which no function of y appears 



g t + [l - D ± /D J E L g x V i g 



(A18) 



At this point it is worthwhile to eliminate between Eqns. (A18) and 

(A15) and obtain 



U i 
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i a/ 
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^i E L 



g„ = (v. + C )g 



a 6 xx (1 - D./D ) B x v i r 
l a 



v . 



or g + 



A 



S„ = TT- 8 



xx \ D - D . / x D 
\ a 1 / a 



(A19) 
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Equation (A19) can be solved with 



g = T(t) e 



ax 



(A20) 



where a is found from 



- e ± Mr + 4 c 2 /D a 



u.e t 

i L 

D-D, 
a i 



(A21) 



Note that D. < D and that £ is anyway a positive parameter. 

1 3 , 

Equation (A21) indicates that the roots are real and there is no periodic 
behavior in x — that striations are not to be found with this description. 

In order for n to oscillate in x , we must have a situation in which 
new constants C^* and C 2 * are introduced such that 



and 



C l* - C 2* = V i 



g t " D a g xx = C l* g 



- D f" = - C * f 
a 2 



(A22) 

(A23) 

(A24) 



Here, Eqn. (A24) will result in (non-periodic) hyperbolic functions in y and 
Eqn. (A23) in oscillatory functions in x provided that 



4C 2 * 




(A25) 



As before we would like to have 



C l* > V i (C 2* = C l* " V 



so as to match L to L_ . 

n T 



(A26) 
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The solution to Eqn. (A24) does not represent the usual diffusion 
controlled situation because in order to meet the symmetry condition at the 
centerline we must give up the physically sensible boundary conditions at the 
wall. Perhaps variable coefficients and the inclusion of recombination can 
remedy this problem with the boundary conditions. Or perhaps we should be 
willing to give up axial symmetry and consider some sort of contraction 
asymmetry which would likely have to include some reconbination because of 
the higher pressures. This, however shall not be pursued in this report. 
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APPENDIX B - INFLUENCE OF TWO-BODY RECOMBINATION OF SCHOTTKY SOLUTION 



The diffusion, equation given in the main text (Eqn. 14) is not suitable 
for elliptic funciton representation, therefore, a different form is developed 
below. 



The diffusion equation with constant coefficients is (without considera- 
tion of whether C^ shall be used) 

.2 2 

- D — y = v . n - an (Bl) 

a , 2 i 

dy 

2 

Defining a stationary state v^n* - an* = 0 , or 

n* = v^/a (b2) 

We can nondimensionalize the equation with 



Obtaining 



n = n/n* and y = y/L 

n 




(B3) 



n' (o) = 0 and n = B , 0 8 1 

n(a/L ) = 0 
n 



It is clear that the boundary condition at the centerline make-up for the 
disappearance of e in Eqn. (B3) . Now Eqn. (B3) can be integrated once to 
yield 



dn 

dy 






2 -3 

3 n 





(B4) 
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or 



(B5) 




The left-hand-side of the above equation can be represented in terms of ellip- 
tic functions. Let 

(n 2 - j n 2 - 8 2 + y B 2 ) = (n - a) (n - b) (n - c) (B6) 

therefore 

(! - e ) W(f - e ) (I + 3s )~ 

a = 2 

b = 8 (0 < S < 1) 

(I - B ) -V(i ' s ) (f + 3S ) 

C = 2 

where a _> b _> n _> c 

With information on a, b, and c we can enter any table for elliptic 
functions (see for example Ref. 18 or 19). Using Byrd and Friedman notation, 
formula 234.00, we have 




where 




dn 

\/(a - h) (b - ft) (h - c) 



= gF 0P,k) 



g = 



2 



V a “ 



c * 



b - c 
a - c 



(B8) 
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and 



'fi = sin 



-1 




(a - c) (b - fl) 
(b - c) (a - n) 



F(<£,k) is tabulated. Note that for 6 = 1 , we get a double root (a = b) 
and the equation is not an elliptic integral (though its solution can still be 
found from the Tables) . When 6=1 we get the boundary-layer behavior 
described by Eqn. (19) or (20) in the main text. 

Except for 6=1, Eqn. B3 represents an eigenvalue problem. In order 
to check this, we solve for the value of the integral at the wall where 
n = 0 . Thus 



As 6 increases, we obtain progressively larger values until at 6=1 we 
get X = «° ) . Figure B1 shows this behavior for the eigenvalue. The results 
shown in Fig. B1 proved identical to results from a computer program at the 
Technion (for the combined diffusion and heat dissipation). 

Increasingly higher degrees of recombination flatten out the density 
profile at the center and increase its slope at the walls. The radial elec- 
tric field is thus enhanced. The variable coefficient solution is not antici- 
pated to change the above results in any major way. 




For very small values of 6 , we obtain the cosine solution, namely 



X = (tt/ 2)^ where L = a/X^^ 

n 
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Figure Bl. Eigenvalue for Increasing Degrees of Recombination 
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Shaped Charges For Neutralization of Buried 
Mines 



A. E. Fuhs, Professor of Mechanical Engineering 
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To develop a technique to neutralize buried 
mines. The technique is to be used by EOD 
officers and crew. 

Mines may be buried at depths from a few inches 
to 6-8 feet. The buried mines are detected and 
located by magnetic means. Error in location 
may be comparable to buried depth. Some means 
of placing an explosive charge near the mine 
is required. A shaped charge can open a large 
hole through which a charge can be rammed. 

An assembly of charges was designed/ built, 
and tested. The assembly of charges consisted 
of a 2-inch diameter shaped charge with % lb 
HE located on top of a cylinder containing 
liquid explosive. The cylinder had a coaxial 
tube which allowed passage of the metal jet 
from the shaped charge. It was found that the 
liquid explosive spray lined the hole; it was 
also found that the cylinder for the liquid 
explosive was rammed down the hole. As a 
result of these early tests, a series of tests 
is being planned using solid explosive located 
below the shaped charge. A time delay deto- 
nator will be used for the solid explosive. 

(See supporting information that follows) 
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NAVAL POSTGRADUATE SCHOOL 
FOUNDATION PROJECT ON UNDERWATER SHAPED CHARGES 

Principal Investigator 
Professor A. E. Fuhs 

Students LT Richard Funk 

LT Larry Kryske 
LT Dana Vugteveen 

As the result of Foundation funding and subsequent additional 
funding from EODF, a series of experiments were conducted at EODF 
Indian Head, MD. The attached material describes the results. 

For FY 78 additional funding was obtained from EODF. It should 
be noted that EODF supports the tests with their own manpower. 

The students who have worked or currently are working on the 
project are as follows: 

LT Richard Funk, USN, LT Larry Kryske, USN and LT Dana 
Vugteveen, USN. 

The attached material is: 

1. Abstract from LT Funk's Thesis. 

2. Invention disclosure on explosive injection by shaped 
charges . 

3. List of contacts at other research and development 
organizations . 

4. Sketch of shock attenuator for mine neutralization weapon. 

5. Description of LT Funk's test at EODF. 

6. Test plan for tests on 22-27 May 1978. 



1 . 

EXPLOSIVE ORDNANCE DISPOSAL MINEHUNTING SYSTEM AND 
MINE NEUTRALIZATION USING SHAPED CHARGE FOLLOW-THROUGH TECHNIQUE(U) 

A brief historical review of mine warfare is given as back- 
ground. The NATO First Scenario serves as the threat scenario 
which provides the impetus for an improved mine countermeasures 
technique. The experience of clearing Haiphong harbor and the 
Suez Canal showed that minesweeping of any sort is a difficult, 
tedious, and lengthy process totally devoid of any glamour. The 
prerequisite of any mine countermeasures operation is training. 

The Navy Explosive Ordnance Disposal community, which is tasked 
with determining the extent of mine fields, has the equipment and 
expertise to initiate this new technique in mine countermeasures. 
This technique entails locating and cataloging all bottom objects 
in United States ports. To counter a mine threat these surveys 
would have to be conducted annually. Experiments were conducted 
which demonstrated the praci ti cal i ty of using the shaped charge 
follow-through concept to neutralize 



buried ordnance. 
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2 . 

INVENTION DISCLOSURE ON 
EXPLOSIVE INJECTION 
BY 



SHAPED CHARGES 



INVENTION DISCLOSURE ON 



EXPLOSIVE INJECTION BY SHAPED CHARGE 

J 

by 

(1 

Richard Funk and Allen Fuha 

There are several military tasks for which the ability to inject an explosiv 
desirable. The military applications are described later. 

Recent tests have shown that either a liquid or solid explosive can be injec 
through sand or mud a considerable distance by use of a shaped charge. The metal 
jet from the charge opens the hole for injection. Figure 1 shows a cross section: 
an arrangement for injecting liquid explosive. Figures 2 and 3 show the apparatu;, 
the shaped charge, in assembled and disassembled condition. Figures 1 and 3 can 't 
compared to identify components. 

Figure 4 is the sketch of apparatus for injecting a solid explosive through ; 

: \ 

sand, or mud. The detonator at the top initiates the detonation wave in the snap, 
charge and starts the time sequence leading to detonation of the high explosive b. 
the shaped charge. This will be discussed in more detail shortly in connection w: 
Figure 6. 

Consider now the sequence of events during the explosion of the shaped chargi 
liquid explosive (LE). At the moment of initiation of detonation, the apparatus : 
sitting on the surface of the sandj see Figure 5(a). The velocity of the metal j< 
shown in Figure 5(b) is about 10000 m/sec. Figure 5(a) shows the condition when 
j®t has bored a hole of 0.2 m or so. The metal jet causes a shock wave in the sal 
Notice the piston which forces the LE from the chamber. Motion of the piston 
relative to the chamber can be ascertained using the reference marks. 

'1 
It 
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In Figure 5(c) the metal jet and slug from the shaped charge have bored a 
large hold in the sand. The chamber is moving into the hole due primarily to the 
hot, high-pressure gas formed by the explosion. A spray of liquid explosive 
develops due to piston motion. In Figure 5(d) the chamber continues to move 
into the hole bored by the slug and the metal jet. The spray of LE coats the sand. 

About a millisecond later the sand begins to collapse as sketched in Figure 5(e). 
The high pressure gases have expanded reducing the pressure on the top of the chamber. 
As a result, the spray starts to go in both directions. Depending on the rate of 
collapse of the sand, the chamber may be trapped part way in the hole. The final 
condition is shown in Figure 5(f). The dark lines are sand soaked with liquid 
explosive (LE). Figure 5(e) is based on observations from tests conducted at EODF, 
Indian Head, MD, July 1977. 

Figure 6 shows a sequence of events similar to Figure 5 but with a solid explosive 
being injected. The injected charge is pushed down the hole by high pressure gases as 
shown in Figures 6(c) and 6(d). By the time the sand collapses, the injected 
charge is at the bottom of the hole as can be seen in Figure 6(e). After a delay 
sufficient to allow complete collapse of the sand, the injected charge is exploded as 
shown in Figure 6(f). 

There are numerous applications which will be summarized below. 

Mine Neutralization . Buried mines can be effectively neutralized. In the 
Navy application, the mines can be buried in mud or sand in a harbor or coastal 
region. In the DSMC application, the mines are of the land variety buried in soil 
jr sand ashore. 

Torpedo Warhead . When used in ASW role, the torpedo would strike a submarine. 

The ballast tank wall would be penetrated followed by the 2 to 3 meters of water 
Inside the tank. The explosive would strike the pressure hull and explode; the 
jressure hull will have been weakened by the metal jet and slug. 
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When used against major combatant surface vessels, the antitorpedo water and 
ftiel tanks along the waterline are penetrated. The solid explosive is injected de 
within the tanks and explodes at that deep location. 

Antipersonnel/Antitank Mine . The apparatus shown in Figure 6(a) could be 
turned upside down and buried 2 to 3 meters in sand, soil, or mud. When detonated 
the metal jet and slug bore a hole upward. Slightly later, the solid explosive 
pops up and explodes. Similar mines exist; however, this mine could be buried ver 
deep (2 to 3 meters). Detection of a plastic, deeply-buried mine would be very 
difficult. 

Anti Bunker . Bunkers frequently consist of a wood or concrete structure with 

0 

sand overlay. The sand may be a few meters in thickness. The shaped charge would 
fired horizontally into the sand; the solid explosive would be injected to a locat 
in contact with wood or concrete. Obviously, the protection offerred by the sand 
negated. 
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Figure 1. Arrangement for injection of liquid explosive. 
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Figure 2. Assembled cylinder for injecting liquid explosive into sand, mud, 
or soil. Hose clamp on left clamps shaped charge. Two 3/8" pipe 
plugs are for filling chamber with liquid explosive. 
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Figure 3. Disassembled cylinder showing components. From left to right across 

the top are the following: two O-rings, disc which shears and becomes 

a piston, tube inserted in end plate, an O-ring, and 3/8” pipe filler plugs. 
From left to right across the bottom are the following: hose clamp, 

hollow cylinder to provide stand-off space for shaped charge, and hollow 
cylinder to hold liquid explosive. 



6 



detonator 




Figure 4. Arrangement for injection of solid explosive. 
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(a) At moment of initiation 
of detonation 



(b) A few 10's of 
microseconds 
later 



(c) A few 100*3 of 
microseconds 
later 



Figure 5. Sequence of events. 



8 



high pressure ga3 




crater 



sand which fell 
in from top 



spray 



liquid 
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pool of 
liquid 
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with 
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(d) A few 100 microseconds ( •) Millisecond (f) Final condition 



Figure 5 (Continued). Sequence of events. 
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high pressure gases 




(a) Initial condition (b) A few 10*3 of (c) Solid explosive 

microseconds entering hole 

later 



Figure 6. Sequence of events for injection of solid explosive. 
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3. 



LIST OF CONTACTS AT OTHER 
RESEARCH AND DEVELOPMENT ORGANIZATIONS 



Individual and Organization 



Extent of Contact 



Colonel Robert Gomez BRL 


Visit 


to 


Laboratory 


Dr. Dietrich 


BRL 


Visit 


to 


Laboratory 


Donald Phillips 


WOL/NSWC 


Visit 


to 


Laboratory 


CAPT Walton Boyer 


WOL/NSWC 


Visit 


to 


Laboratory 


Julius Simmons 


BRL 


Visit 


to 


Laboratory 


William Sykes 


DTNSRDC 


Visit 


to 


Laboratory 


Paul W. Cooper 


Sandia Corp. 


Visit 


to 


Laboratory 


Dr. Salive 


DTNSRDC 


Visit 


to 


Laboratory 


Dr. Julius Enig 


WOL/NSWC 


Visit 


to 


Laboratory 


Jim H'mmond 


NCSL 


Visit 


to 


Laboratory 


Robert Obrachta 


ONR 


Visit 


to 


Office 


Jim Bailey 


ONR 


Visit 


to 


Office 


Dr. V. C. Dawson 


WOL/NSWC 


Visit 


to 


Office 


Frank Romano Code 


034 NAVSEA 


Visit 


to 


Office 



AFAL Eglin AFB 



Visit to 



Colonel 
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SKETCH OF SHOCK ATTENUATOR FOR 
MINE NEUTRALIZATION WEAPON 
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5. 



DESCRIPTION OF LT FUNK'S TEST 



AT EODF FY 77 



DENSITY OF SAND: DRY: 67.2 lbs/cubic foot 

WET: 112.25 lbs/cubic foot 

Tests 1-6 were conducted with the built in stand-off of the 
shaped charge, one cone diameter (c.d.). In all tests other 
than numbers 17-20 a wooden box stand, figure 1, was used to 
support the shaped charge, 

A1 screening was used in tests 1-6 as teller plates the spac- 
ing of the screens was every one inch in the dry then the 
distance the sand settled was noted when the water was added 
and this gave a correction factor for the depth of the screens. 

TEST 1 

Dry sand, 27" in a plastic garbage can, 1 c.d. stand-off 
Penetration: Jet/slug penetrated the bottom of the can, 

leaving a 3/8" diam. hole tnen it penetrated the 1st 3/^" 
oak board of the supporting pallet and cracked the board 
below that one. 

TEST 2 

. Dry sand, 27" in a plastic garbage can, 1 c.d. stand-off 
Penetration: good penetration to the 23" level and then 

burned wood chips found on top of the screen at the 2^" level. 
TEST 3 

Wet sand, depth of sand/water mixture in the can was 25" 
^therefor the correction factor is 0.926. 

Penetration: penetrated the entire can leaving an exit 

approx. 3/8" in diam., hole close to the center of the can 



implying that there was little bending of the jet/slug on 
traveling through the sand/water mixture. 

TEST 4 

Wet sand, settled depth of wet sand 25" therefor the 
correction factor is 0.926. 

Penetration: The jet/slug penetrated the entire can 

leaving an exit hole of approx. lA" and again it was close 
to the center of the can. 

TEST 5 

. Wet sand in a 4 , X4 , X4* plywood box with a correction 
factor of 0,943. A stand-off of 1 c.d. was used. 

Penetration: 28.25" corrected 

TEST 6 

Wet sand in a 4'X4*X4* plywood box with a correction 
factor of 0.943. A 1 c.d. stand-off was used. 

Penetration: 29.25" corrected 

In the tests that follow the number and placement of the 
teller plates was changed. One foot square A1 screens were 
buried at depths of 2 6" and 25" in the dry sand then a one 
foot square 1/8" thick A1 plate (70 75-16) was buried at 24" 
then the cans were filled with sand. Water was then added . 
to the cans which were gently agitated until they stopped 
bubbling. The amount of settling was then noted to be used in 
establishing the correction factor. 
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TEST 7 

2 c.d. stand-off* corrected plate depth 21 . 75 " 

Penetration: Slight dent in the plate about two handfulls 

of burned wood chips resting on the top of the plate. 

TEST 8 

\ 

2 c.d. stand-off, corrected plate depth 21.34” 

Penetration: 3/8" hole in the plate and a pressure rip 

was induced in the bottom of the can. 

TEST 9 

.3 c.d. stand-off, corrected plate depth 22 . 0 ". 

Penetration: severe tearing of the plate resulting in a 

hole 3/4 "X2”, exit hole in the bottom of the can approx. 

3 / 8 ", penetration continued through the 1 st oak board, but 
there was no sign of any penetration of the 2 nd board. 

TEST 10 

3 c.d. stand-off, after noting the depth of original 
settling the can was filled to the top with more sand and 
water this resulted in a corrected plate depth of 24.39". 

Penetration: : 1 / 2 " hole in the plate with very little 
tearing. Bottom of the can ruptured by pressure. _ 

TEST 11 

4 c.d. stand-off, corrected plate depth 21 . 9 " 

Penetration: 1 " hole in plate 

1 / 2 " hole in 1 st oak board 
3 / 8 " hole in 2 nd oak board 
1 / 2 " deep hole in the ground 

TEST 12 

4 c.d. stand-off, added an extension to the can so that 
the distance to the bottom was 30 " instead of 27 ”, corrected 
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plate depth 2^,28". 

Penetration t l/2" hole in plate with some tearing jet/ 
slug missed the pallet boards as it exited the can, but it 
did dent the ground. 

TEST 13 

5 c.d. stand-off, corrected depth to the plate 21.75** 

Penetration: 1 3/4"X3/4" hole in plate 

5/8" hole in 1st oak board 
1/8" hole in 2nd oak board 
1/2" hole in the ground 

TEST 14 

5 c.d. stand-off, refilled can after settling just as in 
test 13, corrected depth to the plate 24.35". 

Penetration: 1/2" hole in the plate caused mostly by 

the tearing action, penetration was to the bottom of the can 
but not through it. 

All remaining test were done in 4 , X4'X4* plywood boxes. All 

\ *T* 

stand-offs totaled 4 c.d. of air before the jet/slug started 
penetrating the wet sand. Teller plates for these tests 
consisted of two A1 plates and seven (7) A1 screen squares. 
These were buried at the original depths of; one plate each 
at 30" and 38" and one square of A1 screening each at 12", 
15*» 18", 21", 24", 27", and 33". The Injector was designed 
at NPS and had a built in stand-off of 3 c.d. through the 
injection chamber. For an exploded veiw of the injector see 
figure 2. 
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TEST 15 

Wet sand with a water head, corrected depth of plate 
27.25" correction factor is 0.909* water depth was 37 H to the 
surface of the sand. S.c. and holder placed in garbage bag. 

Penetration* 3/8" hole in plate, overpressure split 
plate in two, pressure damage only to the screen at 33(30.0") 
TEST 1 6 

Wet sand, corrected depth of plate 27.8" correction 
factor is 0.927. 

Penetration: 1/4" hole in plate with tearing, pressure 

tearing of the screen at 33(30.6"). 

In these last four test the stand-off refers to the distance 
from the base of the cone to the tube running the length of 
the injection chamber. 

TEST 17 

Dry sand, 1 c.d. stand-off (total 4 c.d.), poor technique 
in digging out the box lead to poor results but pictures 
when developed should be consulted. 

Penetration; No red dye found on the surface. At 5.5" 
dye starts and had a circular pattern. Bottom cylinder of 
injector found at a depth of 24" with an 0-ring found at 27" 
and filler plugs found at 30". There was no deformation of 
the plate. 

TEST 18 

Dry sand, 1.25 c.d. stand-off (total of 4.25 c.d) 

Penetration; No red dye on the surface, ring of dye with 
a bullseye presistant to the 18" level then the bullseye 
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starts to move off center, at the 21 M level the ring is still 
6 " in diam. . The bottom cylinder was found at the 24" level 
and sloped down to the 27 " level, the bottom plate off the 
injector was found on the plate at 30 " but it had completly 
turn itself around from what its original position had been. 
The plate was badly cracked, approx. 75/S of its width, and 
the red dye was evedent to the 33 " level. 

TEST 19 

Wet sand with a water head, 1 c.d. stand-off (total of 
4 c.d.). Corrected depth of plate 28.15" and the correction 
factor 0.938. The water depth to the surface of the sand was 
31". 

Penetration 1 The box and tower were scattered all over 
the range so depth of penetration of the different parts was 
impossible to determine. The lower cylinder of the injector 
was recovered. The plate was bowed but other wise unmarked. 
Large amounts of the oil based red dye were found at what had 
been a depth of 42.2" corrected. 

TEST 20 

Wet sand with a water head, 1 c.d. stand-off (total of 
4 c.d.). Corrected depth of plate 27.95" and correction 
factor of 0.932. The water head depth was 36 " to the surface 
of the sand. 

Penetration* Cylinder impacted with the plate and was 
recovered, impact only causing the plate to be dented. Dye 

and small pieces of screen were found on the bottom of the 
box. 
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THINGS TO INVESTIGATE 

1. Use scale model steel cylinders to simulate the mine. 

r 

2. Check penetration of the shaped charge in just water 
using the 4 c.d. air stand-off. 

3. Use different cone material like zinc. 

4. Redesign the injector 

a. make the s.c. holder out of plastic so that the shock 
wave will not be coupled with the base cylinder. 

b. make the base cylinder more aerodynamic ally stable. 

4. Check for optitimum distance between the base of the 
cone and the top of the center tube. 
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EXPLOPED VEIW OF INJECTOR 
FIGURE 2 
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TEST PLAN FOR 

MINE NEUTRALIZATION WEAPON TESTS 



ON 22-27 MAY 1978 
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EUTRALIZATIO 

WEAPON 





TO BE CONDUCTED 
AT 

EODF, INDIAN HEAD 
22.-27 MAY 1978 



Shaped Charge Research Team 
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OF 

MINE NEUTRALIZATION WEAPON 
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SUBMITTED i 

D. VUG TE VEEN LT, USN 

L. KRYSKE LT, USN 
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TEST PLAN FOR EVALUATION OF MINE NEUTRALIZATION WEAPON 
1. Objectives t 

a. Te determine eptimum stand-off fer the Jet Research Center* 
Inc. underwater shaped charge (Y-1076) against sand targets. 

b. T* establish the technical feasibility ef insertion of a 
solid stage explosive into a cavity formed by a shaped 
charge. 



2. Period of Tests > 

21 May 1978 - 27 May 1978 



3« Testing Site * 

U.S. Naval Facility, EOD, Indian Head, Maryland 



4. Synopsis of Schedule of Events : 

20 May 1978 (Saturday) Depart Monterey. PM arrive Dulles; 

proceed to accomodations. 

21 May 1978 (Sunday) Establish contact with EODF Duty 

Officer; liaison with assigned 
contact. Survey test site, install/ 
check-out test equipment. 

22 May 1978 (Monday) Establish contact with -TeUiulUtfT C&- 

Director. Conduct stand-off tests 
against dry sand targets. 



23 May 1978 (Tuesday) 



Conduct evaluation of mine neutralization 
weapon (MNW) against dry sand targets. 



24 May 1978 (Wednesday) 



Conduct evaluation of MNW against 
wet sand targets. 



25 May 1978 (Thursday) 



Conduct follow-up evaluation of MNW 
(wet and/or dry sand targets) 
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26 May 1978 (Friday) Backup day for testing. Prepare 

equipment and samples for return 
to Monterey. 

27 May 1978 (Saturday) Depart for Monterey. 



5. Personnel ! 



Dr. A.E. FUHS, Chairman, Department of Mechanical Engineering, 
Naval Postgraduate School, Advisor 

LT L. KRYSKE, USN Test project officer 

LT D. VUGTEVEEN, USN Test project officer 



6. Requirements : 

a. Test Equipment — see Table 1. 

b. Material — see Table 2. 

c. Support Services 1 

1) Contact personnel for Sunday positioning of equipment. 

2) Photographer f<vr period 22-26 May 1978. 

3) Ordnance personnel for explosive handling. 

4) Hypervelocity range (including use of flash 
radiography equipment). 

7. Safety Precautions ! 

Follow range SOP, no special safety precautions. 

8. Transportation/Handling Hazards ! 

Follow all standard transportation and handling procedures. 



9* Explosive Weights: 

Each test will consist of detonating one one-half pound 
v explosive shaped charge manufactured by Jet Research Center, 

Inc. (Y-1076). 
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10. Disposal Procedures : 

Any excess explosives will be returned to storage. That which 
is in excess and cannot be returned will be destroyed by detonation 
on the final testing day at the discretion of the ordnance safety 
officer. 

11. General Preparations : 

a. Sun 21 May 1978 

1) Locate instrumentation in blockhouse; run all necessary 
cables . 

2) Prepare screens for testing; color code and label. Sizes 
and quantities to prepare: 

3) Attach foil switches to mounting hardware. 

4) Arrange short test rigs for testing: locate on pallets, 

attach eye-bolts and guy-wires; attach pressure trans- 
ducer stand (see drwg ); attach pressure transducer 
housing to first unit (see drwg ); fill all units with 
sand and insert foil switches (see drwg ); attach shaped 
charge positioning hardware (see drwg ). 

5) Time permitting, arrange test apparatus for 23 May testing. 

b. Mon 22 May 1978 (10 shots) 

1) Check out Instrumentation. Position X-ray heads for 
testing (see drwg ). 

2) Place short test rig at hypervelocity site. Attach film 
cassettes to test rig (see drwg ). 

3) Attach instrumentation cables. Insert pressure transducers. 
Position sand bags (see drwg ). 
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4) Insert shaped charge in positioning hardware for stand-off (SO) 

of 1 cone diameter (CD). Insure shaped charge (SC) is 

. 

perpendicular to surface of test rig. Place detonator into 
SC; tape into place. 

5) Steps 2-4 will be repeated 9 more times to determine SC 



penetration at the following SOsi 

2 CD 5 CD 8 CD 

3 CD 6 CD 9 CD 

4 CD 7 CD 10 CD 



Even CD stand-off shots will be conducted with flash 
radiography; odd shots may be fired with or without flash 
radiography depending upon time available for testing. 

6) Prepare long test rigs for testing on 23 May 1978. as before. 




c. 




Tues 23 May 1978 (3 shots) 

1) If necessary continue '•.ests from previous day. 

2) Position X-ray heads for MNW tests, (see drwg ). 

3) Place long test rig at hypervelocity site. Attach film 
cassettes to test rig (see drwg ). 

4) Attach instrumentation cables. Insert pressure transducers. 
Position sand bags (see drwg ). 

5) Insert MNW in mounting hardware. Insure device is perpendicular 
to test rig surface. Place detonator into SC, tape into 
place. 

6) Repeat steps 3-5 for two more shots. 

?) Prepare 3 test rigs for following days' testing using similiar 
procedures as before. Use foil switches designated for 
underwater use. Caulk all openings, fill with water and allow 
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to stand overnight. 

d. Wed 24 May 19?8 ( 3 shots ) 

l) Repeat procedures in steps c - l) - 6) for three shots 



into wet sand targets. 

2) Prepare 3 test rigs for next day's testing (wet and/or dry 
set-ups as determined by test and range officers). 



e. Thur 2 5 May 1978 ( 3 shots ) 

. l) Repeat procedures in steps c - l) - 6) for 3 shots. 

2) Time permitting prepare garbage can rigs for additional 
stand-off data. Fill one can with sand, placing a coded 
teller screen at one inch intervals. Rap can to settle sand 
place last screen on the top of the sand. Invert second can 
(with bottom cut out) and affix to first with ordnance tape. 
Continue sand/screen filling until achieve sand depth of 
4| - 5 feet. 



f . Fri 26 May 1978 ( 5 shots ) 

l) Place SC in positioning hardware as before. Insert detonator 
into SC and tape into place. Repeat for remaining SC. 
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Final Preparations and Test Procedures! 

a. In paragraph 11 all set-ups ended with the placement of detonators 

in the SC. Follow standard hook-up procedures to detonating >-• ■ - 

circuit of hypervelocity range and detonate. For last five tests 
follow standard hook-up procedures to blasting machine and detonate 

b. After range has been certified safe, the test officers will inspect 

all test related items with the range operations supervisor. 

c. After each shot recover as must test apparatus as possible. Remove 
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film cassettes (when applicable) and develop film prior to 
next shot. Determine depth of penetration, record data from 
foil switch electronics. Use polarold cam era to record 
pressure history from storage o ‘scope. Measure hole 
diameters. Recover slug, any pieces of the shaped charge 
and the MNW injector assembly for subsequent analysis at 
KPS. Utilize photo services during this phase, 
d. Due to procedural unknowns and testing environment, the test 
officers are authorized to alter the test setup and pro- 
cedures with the approval of the range safety supervisor. 
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MATERIALS REQUIREMENTS 




MATERIALS REQUIREMENTS (coni) 
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Activation Energies for Elevated 
Temperature Deformation in an Aluminum 
17.5 Weight Percent Copper Alloy 

T. R. McNelley* and A. L. Cipriani** 

ABSTRACT*** 

An Aluminum - 17.5 weight percent copper alloy was prepared 
by casting and subsequent hot rolling to produce a fine dispersion 
of A^Cu in a fine grained aluminum-copper solid solution matrix,. 

This material was evaluated by compressive stress-strain testing 

-4 - 1 

from 300K to 760K, employing strain rates from 1.9 x 10 S to 

- 3 - 1 

9.5 x 10 S . Three regimes in stress were observed for the 

activation energy: between modulus - compensated stresses of 

- 3 -3 

-0.3 x 10 and -0.6 x 10 , the activation continuously increased 

from a value of 107 KJ/mol to 157 KJ/mol ; for modulus compensated 

_ 3 

stress on up to -1.0 x 10 , the activation energy was approxi- 

- 3 

mately constant at 157 KJ/mol; above 1.0 x 10 , the activation 
energy decreased to a constant value of 116 KJ/mol. The strain 
rate sensitivity coefficient in the low stress regime was deter- 
mined to be 0.4, suggesting a tendency toward superplastic deform- 
ation in this material; the observation, however, of the continuous 
varying activation energy in this regime is inconsistent with any 
existing theories for such a deformation process. The activation 
energies at higher stresses, and the stress exponents as well, 
suggest dislocation - climb control of deformation. 

*T. R. McNelley, Assistant Professor of Materials Science, 

Naval Postgraduate School, Monterey, California 93940 

**A. L. Cirpiani, LT, USN, formerly a graduate student. Naval 
Postgraduate School, Monterey, California 93940 

***This paper is being prepared for submission to Acta Metallurgic 
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Warm Working of Eutectic and Eutectoid alloys 
to produce Fine Microstructures 

T. R. McNelley, Assistant Professor of Mechan- 
ical Engineering 

Foundation Research Program (6.2) 

Production of fine, spherodized microstructures 
in Al-Cu and Al-Mg alloys by warm working and 
evaluation of subsequent mechanical properties , 
especially ambient temperature properties. 

The objectives of this research were met. Fine 
two-phase microstructures were produced in both 
an Al-Cu alloy and several Al-Mg alloys. Sub- 
sequent mechanical characterization of these 
alloys has revealed that excellent mechanical 
properties can be developed in Aluminum-Magnes- 
ium alloys containing 11% to 14% magnesium by 
weight. Flow-stress values up to 90,000 psi 
were obtained, for example, in one alloy. This 
has lead to a focusing of this research prima- 
rily on the Al-Mg system; some additional prac- 
tical factors concerning this system are that 
these alloys are less dense than conventional 
aluminum alloys, and the alloys being studied 
also exhibit superplasticity at warm tempera- 
tures. Research will now focus on the process- 
ing and related mechanical properties of these 
Al-Mg alloys. 



A. L. Cipriani, "An Investigation of the Me- 
chanical Properties of Warm-Rolled Aluminum- 
17. 5 weight percent Copper alloy," Master's 
Thesis, December 1976 

F. G. Ness, "High Strength to Weight Aluminum- 
18 Weight percent Magnesium alloy through 
Thermomechanical Processing," Master’s Thesis, 
December 1976 

T. L. Glover, "Effects of Thermo Mechanical 
Processing on Aluminum-Magnesium, Alloys Con- 
taining High Weight Percentage Magnesium" 
Master's Thesis, December 1977 



C. P. Bingay, "Microstructural Response of 
Aluminum-Magnesium Alloys to Thermomechanical 
Processing," Masters Thesis, December 1977 
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ABSTRACT 

Mechanisms of microstructural refinement in aluminum- 
magnesium alloys were investigated. Alloys containing from 
15% to 19% magnesium were exposed to various processing 
schemes, and the resulting microstructures were examined. 
Isothermal forging resulted in some refinement depending on the 
temperature, strain and strain rate. However, all samples iso- 
thermally forged exhibited microstructures having relatively 
large amounts of the brittle intermetallic phase present. 

The addition of tertiary alloy elements resulted in little im- 
provement in microstructure. Non-isothermal schemes resulted 
in the most promising microstructures. A process of high 
temperature soaking followed by deformation at relatively 
lower temperatures is concluded to be the most promising 
means for processing alloys with magnesium contents below 
15% by weight. 



Bingay, Charles P., "Microstructural Response of Aluminum- 
Magnesium Alloys to Thermomechanical Processing," Master's 
Thesis, December 1977. 



ABSTRACT 



An aluminum-17. 5 weight percent copper alloy was 
warm rolled to achieve refinement of the 
oicrostructur e. This refined microstructure consisted 
of finely dispersed intermetallic A^Cu particles with 
an average size of 1.2 microns in an aluminum matrix. 
This led to improved room temperature properties as 
well as the onset of superplasticity at elevated 
temperatures. Ductility and toughness were increased 
almost six-fold, despite slightly decreased yield 
strength and maximum compressive strength. It would 
appear that warm working can be used to enhance both 
room temperature properties and superplasticity, and 
that further research could maximize this improvement. 



Cipriani, Alfred L., "An Investigation of the Mechanical 
Properties of Warm Rolled Aluminum-17.5 Weight Percent 
Copper Alloy," Master's Thesis, December 1976. 
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ABSTRACT 

Cast ingots of three Al-Mg alloys, containing 7% Mg, 

151 Mg, and 19% Mg, were obtained from Kaiser Aluminum and 
Chemical Corporation Center for Technology. Billets were 
machined from these castings and upset forged at . 9 % of 
either the solvus or eutectic temperature, as appropriate, 
to a true strain of 1.5. These alloys were then evaluated 
by compressive s tress - strain testing at various temperatures 
and strain rates to determine both the ambient and elevated 
temperature characteristics. Due to the inability to produce 
very fine second phase particles by upset forging, the 15% 
and 19% Mg alloys were high in strength but brittle at room 
temperature. However, refinement of the second phase parti- 
cles during compression testing led to superplastic behavior 
at elevated temperature, with a strain rate sensitivity 
coefficient of .43 being achieved in the 19% Mg alloy. It 
was further observed that the addition of Mg increased strain 
rate sensitivity in these Al-Mg alloys at all temperatures, 
especially at elevated temperatures. 

Glover, T. L. , "Effects of Thermo-Mechanical Processing on 
Aluminum-Magnesium Alloys Containing High Weight Percentage 
Magnesium," Master's Thesis, December 1977. 
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ABSTRACT 



An Aluminum- 18 weight percent Magnesium alloy was 

o 

prepared by casting and then warm rolling at 425 C to 
94% true strain. This alloy was compression tested at 
six strain rates from 0.00664 per minute to 0.332 
per minute and at eight temperatures ranging from 25* C 
to 42 ^ C. The most significant result is that a warm 
rolled Aluminum-18 weight percent Magnesium alloy can 
exhibit compressive strengths in excess of 95 ksi, in 
a material of 10% lower density than commercial high 
strength Aluminum alloys. Furthermore, one can 
envision a thermal mechanical process involving warm 
working followed by cold working at room temperature 
whereby one can attain an ultimate tensile strength 
greater than 90 ksi. Additionally, superplastic 
behavior at elevated temperatures was manifest in the 
relatively high value of strain rate sensitivity and 
by the value of the activation energy for deformation. 



Ness, Frank G. Jr., "High Strength To Weight Aluminum- 
18 Weight Percent Magnesium Alloy Through Thermal Mechanical 
Processing," Master's Thesis, December 1976. 
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Interval Modulation of a Sinusoidal Carrier 

G. Myers, Associate Professor of Electrical 
Engineering 

Foundation Research Program (6.2) 

This work is concerned with radio commu- 
nications. This research investigated use 
of a modulated interval, between bursts of 
a sinusoidal carrier, to convey the infor- 
mation in a message waveform. We called 
this form of modulation “interval modula- 
tion" (IM) . 

IM produces a delay between bursts of the 
sinusoidal carrier. This delay is propor- 
tional to the amplitude of the modulating 
message waveform. Thus, IM is unlike AM, 

PM, or PM where the amplitude, frequency 
or phase of a sine wave vary in accordance 
with a message waveform. 

The research resulted in mathematical des- 
criptions of the IM carrier in the time 
domain and frequency domain. Early in the 
investigation, we completed the design of 
an IM modulator and demodulator. We then 
built and tested these circuits. The tech- 
nique was then verified by transmitting 
voice and music signals through the system. 
The quality of the recovered signals was 
excellent. 



G. Myers, "Interval Modulation (IM) of a 
Sinusoidal Carrier", Technical Conference, 
SIEEM 77, Monterey, Mexico, September 77. 

G. Myers, "Symposium on Spread Spectrum 
Communications", Naval Postgraduate School, 
September 1976. 
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G. Myers and E. L. Kilborn, Jr., "Interval 
Modulation (IM) of a Sinusoidal Carrier", 
Technical Report, NPS62Mv77031, February 
1977. 

E. L. Kilborn, Jr. , "Interval Modulation 
(IM) of a Sinusoidal Carrier", Master's 
Thesis, December 1976. 
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ABSTRACT 



INTERVAL MODULATION (IM) CARRIER SIGNAL 

Glen A. Myers Edgar L. Kilborn, Jr. 

Naval Postgraduate School Naval Postgraduate School 

Code 6 2Mv SMC Box #1657 

Monterey, CA 93940 Monterey, CA 93940 

In classical radio communications, either the amplitude A, 
the phase 4>(t), or the frequency - of a carrier wave A cos 

<j>(t) is caused to vary in accordance with the message m(t) 
transmitted. This work uses a time delay to modulate a carrier 

; 

in accordance with m(t) rather than amplitude (AM) or phase ( Pb 
or frequency (FM) modulation of the carrier. 

Operation of the time delay or interval modulation (IM) if 1 
as follows: the message m(t) is sampled and after a delay of t 
seconds, where t is proportional to that sample value of m(t), 
a burst of sinusoidal carrier is transmitted until the next sai 
pie of m(t) is taken. The sequence then repeats to define the 
IM carrier signal. Each burst of the carrier signal is of fix* 
duration (constant number of cycles transmitted per burst) theii 
by making the sampling of m(t) non-periodic. 

The IM carrier signal is useful only if the message (sampjf 
can be recovered from the intervals t between bursts in a rece 
Present work involves the use of a phase-locked loop (PLL) as 
IM carrier signal demodulator. To date, successful generation 
and detection of IM carrier signals has been accomplished usin 
voice and music as the modulating (message) signal. 
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[ NTERVAL M ODULATION (IM) OF SINUSOIDAL CARRIER WAVES 

by 



Glen A. Myers 

U.S. Naval Postgraduate School 
Monterey, California 

In classical radio communications, either the amplitude A, the phase 
4>(t), or the frequency of a carrier wave A cos <f>(t) is caused to 

vary in accordance with the message m(t) transmitted. In this paper we 
consider use of a time delay between bursts of a pure sine wave to carry 
the message m(t). This technique, which is referred to as interval modu- 
lation (IM), is distinct from amplitude (AM), or phase (PM) or frequency 
(FM) carrier modulation. 

The principle of IM is explained using waveforms, and parameters 
of the IM carrier are defined. Then a frequency description is considered, 
and the relation of modulated carrier bandwidth to system parameters is 
discussed. 

A block diagram of the modulator used in the experimental work is 
provided, and its operation is briefly treated. A conventional phase- 
locked loop (PLL) is used to demodulate the IM carrier. Voice and music 
used to modulate the carrier were recovered using the PLL. The quality 
of the demodulated voice and music is very good. 

Applications are briefly mentioned. 

1 . Time Waveforms 

Interval Modulation (IM) of a sinusoidal carrier conveys the informa- 
tion of a modulating message waveform m(t) by sampling m(t) at the comple- 
tion of burst of the sinusoidal carrier and varying the interval before the 
next burs!- oy amount proportional to the amplitude of the sampled value 
of m(t; 1 r ■ /en oc voltage used to modulate the carrier, the intervals 



between the bursts of the sinusoidal carrier are fixed ano equal. There- 
fore, the IM signal modulated by a constant voltage i\ pei iodic. An 
example of the IM signal for a particular DC modulating voltage is shown 
in Fig. 1 where is the duration of the burst and is the interval 
between bursts. 

For a time-varying (AC) modulating voltage m(t), the interval 
between bursts of the sinusoidal carrier is caused to vary according to the 
amplitude A of the sample of this modulation or message voltage m(t). 

In all the experimental work, T is caused to vary linearly with A. 

A 

IM is closely related to pulse width modulation (PWM), which is also 
referred to as pulse duration modulation (PDM). Let Vp^(t) represent PWM 
as shown in Fig. 2. Then, we can write the IM carrier which we call Vj^(t) 
as 

v IM (t) = [ 1 - v pw (t)][V sin 2v f c 1 3 

where V and f are the amplitude and frequency of the sinusoidal carrier. 

2. Frequency Description 

The Fourier transform of Vj^(t) is the frequency description of inter- 
est. The equation for Vj^(t) previously written in terms of Vp^(t) shows 
that the spectrum of the IM carrier consists of a component at the carrier 
frequency and sidebands which are the frequency description of Vp^(t), the 
PWM signal . 

For an arbitrary message m(t), the Fourier transform of Vp^(t) is 
not available. Results are available for particular periodic voltages m(t). 
For example, the Fourier transform of Vp W (t) when m(t) - M cos 2uf m t is 
presented in Black.* For even this simple case, however, the transform 

H. S. Black, Modulation Theory , D. Van Nostrand Company, Inc., 1953, p. 275, 
Eqs . 17-25. 
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is a complicated function containing many side frequencies separated by 

f s Hz ,,nd havin 9 magnitudes given by the ordinary Bessel function. Here, 
f $ is the sampling rate. 

Laboratory work confirms the presence of many side frequencies. Fig. 3 
is the spectrum of the voltage shown in Fig. 1. in this case, m(t) is a 
constant (DC voltage) . 

For an arbitrary m(t), the side frequencies caused by the sampling 
are present. Each of these side frequencies as well as the carrier appears 
to be modulated by m(t). The local spectrum about the carrier and each side 
frequency is similar to that generated by FM. That is, the message m(t) 
appears to individually frequency modulate each spectral line. 

An example of the spectrum of the IM carrier when m(t) is a tone 
(sine wave) of frequency f m Hz is shown in Fig. 4. In this case, the 
carrier and four side frequencies are shown. The amplitude of m(t) is 

large, and this results in many lines separated by f m Hz surrounding the 
carrier and each side frequency. 

Tim nature of the spectrum of the IM carrier is of interest 
for at least two reasons. First, the broadband character may have applica- 
tions as mentioned later. Second, the presence of an FM type of spefctrum 
about each side frequency suggests a method of demodulating the IM carrier. 

3. The Modulator 

The IM signal modulator uses samples of the input message waveform 

m(t) to determine the length of the interval between the bursts of the 
sinusoidal carrier. 

'he for ft. of the IM signal depends upon several parameters. These 
parameters include the number of cycles in a burst, the frequency of the 
sinusoidal racier comprising the burst, the length of the basic interval 
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between bursts, and the range of variation (due to the amplitude of the 
modulating signal) of the interval between bursts. For this reason, the 
IM signal modulator is designed to permit, within certain limits, varia- 
tion in the IM signal parameters. 

Fig. 5 shows the basic block diagram for the IM signal modulator. 
Referring to Fig. 5, the frequency of the sinusoidal carrier is varied by 
changing the frequency of the Gated Sinusoidal Burst Generator. The number 
of cycles of the sinusoidal carrier in each burst is controlled by chang- 
ing the length of the Burst Window Generator. The basic interval between 
bursts is controlled by changing the slope of the ramp output from the 
Gated Ramp Generator. And, the range of variation of the interval, due 

to modulation, is controlled by the amplitude of the modulator's input 
signal . 

Fig. 6 is a simplified timing diagram for the IM signal modulator and 
shows the sequence of events which generate the IM signal. The operation 
of the modulator is as follows. Coincident with the completion of a burst 
of the sinusoidal carrier a sample command is issued to the Sample and Hold 
circuit and the Gated Ramp Generator is turned on. The outputs from the 
Sample and Hold circuit and the Gated Ramp Generator are fed through a 
Differential Amplifier having zero volts output when the output of the 
Sample and Hold circuit is equal to the output of the Gated Ramp Generator. 
The Zero Crossing Detector following the Differential Amplifier notes the 
zero crossing and turns on the Burst Window Generator. The Burst Window 
Generator then gates on the Gated Sinusoidal Burst Generator. The length 
of the burst from the Gated Sinusoidal Burst Generator is a function of 
the length of the Interval Modulated Gate signal produced by the Burst 
Window Generator. The length of the Interval Modulated Gate signal is 
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controi led by a tuning sequence, the duration of which is variable. Upon 
completion of the timing sequence the Interval Modulated Gate signal is 
turned oft, another sample is taken, and the sequence of events repeats. 

4. The Demodulator 

The time representation of the IM carrier indicates that samples of 
the message can be recovered with a circuit that converts the duration of 
the interval between bursts to a corresponding voltage. This task was of 
concern early in the work because of the microsecond variations expected 
for audio messages m(t). It was decided to investigate use of a phase- 
locked loop (PLL) as a possible demodulator. The notion was that the 
transient behavior of the PLL between bursts of the carrier might provide 
demodulation. 

The PLL turned out to be a fortunate choice. A standard PLL chip 
recovered the message from the IM carrier. Further, it was discovered 
that the PLL could be "tuned" to many different discrete frequencies in 
the region about the carrier and recover the message. 

The nature of the spectrum of the IM carrier indicates why a PLL can 
be used as a demodulator. Further, the reason demodulation is possible 
at many frequencies is explained by the presence of discrete side frequen- 
cies. In this application, then, the frequency description rather than the 
time description was found to be more useful in considering the demodula- 
tion process. 

Limitations on the use of the PLL as a demodulator arise when the 
bandwidth of the modulated carrier (or a side frequency) exceeds the spac- 
ing between frequencies (sampling rate). This occurs when the ampli- 

tude of ffie > ?•••!,; • ing signal is large. (This corresponds to large modula- 
tion iruK- • 11 ■ pm.': To prevent such overmodulation, reduce the level of 

the mod ; : < • . > . 
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4 . Application s 

The work presented here was not motivated by any particular problem 
nor dictated by any specific application. The nature of the spectrum of the 
IM carrier suggests some possible uses. 

The IM technique generates several "sub-carriers" of controlled 
spacing. Each of these sub-carriers is modulated (angle modulation) 
by the message m(t). This means a single transmitter can be used with 
several receivers each tuned to a different frequency. 

Alternatively, any one receiver can select the sub-carrier which 
provides the best reception. The IM technique may be useful, then, where 
use of frequency diversity is desirable. 
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Figure 1. IM SIGNAL MODULATED BY A CONSTANT VOLTAGE. 

The carrier frequency is 200 kHz, T = lOysec, 

and ^ = 20 .sec . x 



Figure 3. 



I inure 4. 
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SPECTRUM OF IM SIGNAL MODULATED BY A CONSTANT 
VOLTAGE. The frequency range of 0 to 450 kHz 

;s shewn. 




SPECTRUM OF THE IM 
MODULATING SIGNAL, 
is a sine wa ve . 



SIGNAL for LARGE AMPLITUDE 
The modulating signal m(t) 
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Computer Aided Acoustical Imaging 
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Summary : 



J. P. Powers, Professor of Electrical Engi- 
neering 

Foundation Research Program (6.2) 

This work is part of a continuing effort to 
explore the use of computer processing of co- 
herent (both amplitude and phase information 
present) data to generate images of objects 
from transmitted or reflected ultrasonic 
fields. The short term objective was to con- 
struct an experimental system to record the 
complex data fields of an ultrasonic wave for 
computer processing. Additionally it was 
desired to explore the presentation of the 
data on a computer driven video display for 
presentation of the data or resulting image. 

The first acoustic images using this system have 
been obtained. The data aquisition system has 
been built and successfully tested. Prelimi- 
nary acoustical images have been displayed on 
a Ramtek Gx-100 television display with 16 
gray levels (limited by the display refresh 
memory capability) . The data acquisition sys- 
tem consists of a logarithmic compressive am- 
plifier with 60db of dynamic range, a phase 
detector, and a high accuracy positioning and 
recording system. Testing and calibration of 
the hardware implementation of these devices 
shows acceptably high degrees of accuracy and 
repeatability. The data is recorded on analog 
tape. After digitization of the recorded data 
and subsequent demultiplexing of the channels 
the data has been successfully stored in mem- 
ory of PDP 11/40 computer for subsequent dis- 
play or processing. Using the Ramtek Gx-100 
display driven by a PDP 11/40, interactive 
display capability has been developed to allow 
the operator to call up displays of amplitude, 
phase or intensity data. He also has the capa- 
bility to change the quantization levels, the 
gray scale assignment and to assign various 
pseudocolor values to the data. This inter- 
active capability was designed and tested dur- 
ing the research period using both computer 
generated and experimentally recorded data. 
Future efforts include further expansion of this 
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Publications : 



interactive display capability, refinements 
in the data acquisition system and investiga- 
tion of more experimental objects to assess 
the overall capabilities of the system. 

J. P. Powers, "Computer Simulation of Linear 
Acoustic Diffraction, Acoustical Holography ,” 
Vol. 7, L. W. Kessler, Ed., Plenum Press, New 
York, 1977, pp. 193-205. 

J.P. Powers, Lt. R. T. O'Bryon, Lt. J. W. 
Patton, "Ultrasonic Imaging System Incorporat 
ing Computer Aided Coherent Processing", Ab- 
stracts, Ninety-fourth Meeting of the Acousti 
cal Society of America, Miami Beach, 12-16 
December 1977, in J. Acoustical Society of 
America , 62 (Supplement 1): S22, 1977. 



COMPUTER SIMULATION OF LINEAR ACOUSTIC DIFFRACTION 



John P. Powers 

Department of Electrical Engineering 
Naval Postgraduate School 
Monterey, California 93940 



ABSTRACT 

Computer-aided acoustical imaging systems and computer simula- 
tions of other acoustic imaging techniques frequently require simu- 
lation of linear acoustic diffraction of large complex- valued data 
arrays. Computation efficiency requires the use of fast Fourier 
transform techniques. This paper compares two Fourier transform 
formulations of the propagation problem: the Fresnel integral and 

the spatial frequency domain approach. The following features are 
compared: restrictions on maximum and minimum propagation distances, 

sample sites and number of samples required, adaptability to image 
processing techniques, and computational requirements. 



INTRODUCTION 

The use of computers in computer-aided acoustic imaging has be- 
come increasingly popular in recent years. The use of the computer 
in obtaining images by such techniques as backward wave propagation^ 
offers such advantages as the elimination of the reconstruction wave- 
length scaling problem^ obtained with optical reconstruction tech- 
niques, reference- free holography that uses the linear detection 
properties of piezoelectric transducers, and the possibility of in- 
corporating image enhancement to improve the image obtained. Addi- 
tionally computer simulation of the holographic process has been a 
useful tool in studying such novel techniques as phase-only holo- 
grams' 5 and kinoforms 4 . In most computer-aided imaging techniques it 
is necessary to mathematically simulate the scalar wave diffraction 
process. This paper compares two techniques for this simulation. 



describing their features and some of their advantages and disadvan- 
tages . 

In representing the acoustic diffraction formulations we will 
implicitly assume that the propagation medium is linear and homoge- 
neous. We also will require that the resulting diffraction inte- 
grals must be amenable to computer solution in a reasonable amount 
of time which at the present infers that the integrals must be in 
the form of Fourier transforms so that the speed and efficiency of 
the Fast Fourier Transform^ (FFT) can be brought to bear on the 
problem. The general problem then is: given a complex scalar wave 

U^(xi, yi, 0 ) at some input plane, find an expression for the wave 
Uq(xo, Yo> 2 ) at some parallel output plane a distance z away, subject 
to the wave equations. Two forms of the solution incorporate the 
Fourier transform and will be considered after a short review of the 
features of the analog Fourier transform and the discrete Fourier 
transform (DFT) . 

The two dimensional analog Fourier transform is defined by the 
relationships 

U(x,y) = f2 A(u,v)e^ 2 ^ UX+v> - ) dudv =^*' 1 {a(u,v)[ (1) 

A(u, v) = // U(x,y)e" j2rCuX+vy) dxdy =^|u(x,y)} (2) 

where LJ(x,y) is the complex function in the space domain; A(u,v) 
is the complex transform in the spatial frequency domain; 
u,v are spatial frequencies (dimensions of cycles/meter) ; and 
3“, 3" -1 are symbolic operators for the transform and in- 
verse transform operations respectively. 

The discrete version of the Fourier transform (assuming an equal 
number of samples and sample spacing in both dimensions) is given by 

f(ma,na) . 1 £ £ ACkK,in)e j2, 'C Cl ’ a '“ ,+(na - £2 fl 

N^k=o £=o 

= F~*j a (kn,£n)[ (3) 

A(k Q,£fl) = J]* ^ f(ma,na)e"^ 2ir 

m=o n=o 



= F 



| f (ma,na) } 



(4) 



where o<_ m,n£ N-l; 
o< k,£< N-l; 



